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ABSTRACT 


The design of a high resolution radar for altimetry and ocean 
wave height estimation is studied. From basic principles, it is shown 
that a short pulse wide beam radar is the most appropriate and recommended 
technique for measuring both altitude and ocean wave height. To achieve 
a topographic resolution of + 10 cm RMS at 5.0 meter RMS wave heights, as 
required for SEASAT-A, it is recommended that the altimeter design include 
an onboard adaptive processor. The resulting design, which assumes a 
Maximum Likelihood Estimation (MLE) processor is shown to satisfy all per- 
formance requirements. A design summary is given for the recommended radar 
altimeter, which includes a full deramp STRETCH pulse compression technique 
followed by an analog filter bank to separate range returns as well as the 
assumed MLE processor. 

A feedback loop implementation of the MLE on a digital computer 
is examined in detail. Computer "size 11 , estimation accuracies, and bias 
due to range sidelobes are given for the MLE with typical SEASAT-A param- 
eters. The standard deviation of the altitude estimate is developed and 
evaluated for several adaptive and nonadaptive split-gate trackers. Split- 
gate tracker biases due to range sidelobes and transmitter noise are exam- 
ined. An approximate closed form solution for the altimeter power return 
is derived and evaluated. 

The feasibility of utilizing the basic radar altimeter design 
for the measurement of ocean wave spectra as well, is examined. A prelim- 
inary design analysis shows the resulting system parameters for the ocean 
wave spectrometer to correspond very closely to those of the radar altim- 
eter design. The required modifications are that the transmitter output 
is switched to a 2 meter antenna which is steered 20° away from nadir and 
scanned conically with a period of 5 seconds. On receive, the full deramp 
processor is replaced by a surface wave pulse compressor and detector fol- 
lowed by a bank of 11 one-third (1/3) octave filters to perform the spectral 
analysis. The filter outputs are then detected and integrated for a period 
of 26 fisec per pulse. This preliminary design meets all specified perform- 
ance requirements . 
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Part I 

ALTIMETRY AND OCEAN WAVE HEIGHT 

1.0 INTRODUCTION AND SUMMARY 

The design of a high resolution radar for altimetry and ocean 
wave height estimation is examined in this part of the report. To achieve 
a topographic resolution of 10 cm RMS over practically all sea surface 
roughness, as required for SEASAT-A, it is recommended that the design include 
an on-board adaptive processor. The resulting design, which assumes a 
Maximum Likelihood Estimation (MLE) processor, satisfies all the specified 
SEASAT-A performance requirements for altitude and wave height measurements. 

Section 2.0 examines the basic principles involved in the measure- 
ment of altitude and wave height. Consideration of both the short pulse 
radar and two frequency interferometer technique with either a wide or nar- 
row antenna beam, leads to the conclusion that a short pulse wide beam radar 
is the most appropriate and recommended technique for measuring both altitude 
and wave height. A design summary is given for the recommended radar which 
includes a full deramp STRETCH pulse compression technique followed by an 
analog filter bank to separate range returns as well as the assumed MLE 
processor. 

Section 3.0 describes the implementation and performance of a 
joint MLE which would simultaneously estimate altitude, wave height, and 
signal-to-noise ratio. A step-by-step description of a feedback loop imple- 
mentation of the MLE on a digital computer is presented to show the simplic- 
ity of the concepts involved. It is shown that a relatively modest mini- 
computer would be of adequate "size" for implementing the processor. Esti- 
mation accuracies as well as biases associated with range sidelobes are 
given for typical SEASAT-A parameters. Performance requirements are shown 
to be achieved with a signal-to-noise ratio greater than 5 dB even at 20 
meter significant wave height. The biases due to range sidelobes are shown 
to be negligible when the average RMS sidelobes are less than -50 dB. This 
sidelobe level can be achieved in practice. 
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Section 4.0 examines the performance of a split-gate tracker. An 
easily interpreted expression for the variance of a split-gate tracker is 
derived and used to show the effect of changing gate widths and placement 
of the early gate. A comparison of tracker accuracies for several adaptive 
and nonadaptive split-gate trackers is given for typical SEASAT-A parameters. 
Altitude bias due to range sidelobes is examined. It is shown that while 
receiver weighting to reduce peak sidelobes ’will not be necessary, the aver- 
age far out sidelobe level must be less than -50 dB RMS in order to keep the 
altitude bias less than 6 cm. It is also shown that an altitude bias is 
caused by the additive transmit noise burst resulting from gating the TWT 
on prior to (and keeping in on after) transmitting the linear FM signal. 

To maintain the bias due to this effect at an acceptable level, the ratio 
of total transmit signal energy to total transmit noise energy must be at 
least 17 dB. 

Section 5.0 derives an approximate closed form solution for the 
shape of the altimeter mean power return. The radar equation for a distrib- 
uted target is derived and applied to the particular case of a satellite 
altimeter. The resulting closed form solution is given in terms of the 
usual radar parameters, an "effective" range distribution, which is a com- 
bination of RMS wave height and RMS range resolution, and an "effective" 
beamwidth which is a combination of ocean surface slope distribution and two- 
way antenna beamwidth. Illustrations showing the shape of the mean power 
return as a function of RMS wave height and "effective" beamwidth are given 
for typical SEASAT-A parameters. It is shown that a small "effective" beam- 
width causes a distortion in the shape of leading edge of the mean return and 
produces a wave height dependent altitude bias. To eliminate this bias, an 
antenna beamwidth of 5° to 10° would be required for SEASAT-A. Such a small 
antenna would, however, require a corresponding increase in peak transmitter 
power and/or compression ratio to achieve the recommended signal-to-noise 
ratio. It is recommended that a study be made of all of the trade-offs 
involved (including the effects of antenna pointing error) in order to deter- 
mine an optimum antenna beamwidth. 
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Appendix A presents the theoretical derivations for a broad class 
of estimators which are obtained by minimizing a penalty functional that 
depends on both the observed data and on the parameters to be estimated. 
Maximum Likelihood Estimators (MLE) are special cases of this class of esti- 
mators. Formulas for computing the asymptotic bias and variance of the 
estimators are derived. A feedback loop implementation scheme is described 
and formulas for computing the mean response and the variance of the loop 
estimates are derived. The general solutions are then evaluated for the 
satellite altimetry case. The resulting formulas provide the basis for the 
results described in Section 3.0 of the report. 
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2.0 


HIGH RESOLUTION RADAR IMPLEMENTATION 


It is desirable, for geodetic purposes as well as the measurement 
of many oceanographic phenomena of interest, to measure the position of the 
mean sea surface to an accuracy of +10 cm RMS. The design of a satellite 
radar altimeter with this accuracy would be quite simple were it not for the 
unavoidable degradation in performance caused by sea surface roughness. 

Recent studies ^ have shown that the standard deviation of the altitude 
measurement is directly proportional to sea surface roughness. Thus an 
altimeter design having a 10 cm RMS accuracy at 2 meter wave heights (peak- 
to- trough) would only provide a 100 cm RMS accuracy at 20 meter wave height. 

The design of a satellite altimeter having a topographic resolution of 10 cm 
RMS over practically all sea surface roughness (1 to 20 meters peak-to-trough) 
cannot be achieved on board the spacecraft with standard "split gate" tracking 
algorithms . 

To circumvent this problem, it is recommended that the altimeter 
design include an on board adaptive processor. At the present time several 
candidate adaptive processors are under consideration. These include a 
Maximum Likelihood Estimate (MLE) processor , a Minimum Mean Square Error 
(MMSE) processor and several Adaptive Split-Gate (ASG) trackers 
which either continuously or discretely approximate the weighting functions 
of the MLE or MMSE processors. The theoretical basis for both the MLE and 
MMSE processor was developed by Technology Service Corporation for NASA/Wallops 
under Contract No. NAS6-2241. From a performance point of view, the MLE is best 
since it would provide optimum (minimum variance) estimates of satellite 
altitude, ocean wave height, and electromagnetic ocean surface reflectivity. 

As such, only the salient parameter estimates need to be transmitted to the 
ground instead of the large amount of raw data otherwise required for ground 
processing. The MMSE processor and continuous ASG tracker represent pro- 
gressive simplifications, respectively, in processor design and complexity 
with a corresponding reduction in performance. Finally, the discrete ASG 

tracker approximation to the weighting functions of the MMSE processor proposed 

( 2 ) 

by MacArthur represents the simplest processor design and complexity with 
a corresponding performance that marginally meets the requirements for SEASAT-A, 
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In addition to the obvious trade offs between performance and complexity, 
the recommendation of a particular on board adaptive processor for SEASAT-A 
requires further study and evaluation of the many subtleties associated with 
each processor. However, in order to specify a complete system design, the 
remainder of this section of the report assumes the selection of an on board 
MLE processor. 

Section 2.1 examines the basic principles envolved in the measure- 
ment of altitude and wave height. Consideration is given to both the short 
pulse radar and two frequency interferometer techniques operating with either 
a wide or narrow antenna beam. The ensuing discussion leads to the conclusion 
that a short pulse wide beam radar is the most appropriate and recommended 
technique for measuring both altitude and wave height. 

Section 2.2 summarizes the design of the recommended radar which 
includes an on board MLE processor. The recommended pulse compression 
technique, because of bandwidth considerations, is a full deramp STRETCH 
followed by an analog filter bank to separate range returns. The resulting 
design achieves all specified performance requirements for altitude and wave 
height measurements. 

2 . 1 Basic Principles and Recommended Approach 

The basic purpose of a satellite altimeter is to measure the height 
from the satellite to the mean sea surface. For a nadir looking altimeter 
this requires the measurement of the mean of the density of specular point 
scatterers with respect to height. A measurement of wave height may also 
be obtained with the altimeter if the spread (standard deviation) of the 
density of specular point scatterers with height is available. Rigorous 
theoretical developments of techniques for achieving altitude and wave height 
measurements are readily available in the literature^ . For the most 
part these works by the very nature of their rigorous and precise developments 
tend to obscure the simplicity of the concepts involved. The following is an 
attempt to explain the principles envolved at the expense of rigorous 
justification. 
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First consider a short pulse radar altimeter which, depending upon 
the geometry, can be operated in either a narrow antenna beam or wide antenna 
beam mode. The impulse response (mean power response of a matched filter 
receiver) of a smooth flat ocean would correspond to the radar cross section 
versus range, <7(R), and would appear as shown in Figure 2.1-b. In the wide 
beam (pulse limited) mode cr(R) is a step response which occurs at the location 
of the surface and lasts until the intersection of pulse and surface passes 
out of the antenna beamwidth. In the narrow beam (beam limited) mode, the 
c? (R) is simply the transmit pulse with an epoch corresponding to the location 
of the surface. Next consider a discrete distribution of specular points with 
height as shown in Figure 2.1-c. Now the composite return signals are simply 
the superposition of impulse responses weighted by the distribution function. 

Thus, in the wide beam mode the leading edge of the impulse response 
corresponds to the integral of the specular point density (specular point 
cumulative distribution function). Here an altitude measurement can be obtained 
by tracking the half power point of the leading edge of the return signal (the 
median) and an estimate of wave height is obtained by measuring the rise time 
(shape) of the leading edge. In the narrow beam mode the impulse response 
corresponds to the specular point density itself. Thus, altitude measurement 
is obtained by tracking the centroid of the return signal and an estimate of 
wave height is obtained from the width of the return signal. Passing now to a 
realistic model for the specular point density of the sea surface, Barrick^ 
has shown to a first order approximation that the density is normal. Thus the 
impulse responses for the wide and narrow beam modes would be as shown in 
Figure 2.1-d. 

Now a comparison of the two modes of operation at typical satellite 
altitudes shows the narrow beam mode to be impractical because of the large antenna 
required. That is, in order to have an impulse response which adequately 
corresponds to the wave height density, the differential path length at the 
edge of the beam &R, is required to be less than the smallest wave height of 
interest and the range resolution must be about half the smallest wave height 
of interest. For a 1 meter wave height and satellite altitude of 725 km this 
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Figure 2.1. Impulse Response of a Short-Pulse Radar For Narrow- 
and Wide -Beam Antenna. 
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translates to an antenna diameter of at least 20 meters and a pulse width of 
about 3 nsec. On the other hand, the wide beam mode simply requires that the 
differential path length AR be much greater than the largest wave heights of 
interest and the range resolution be about half the smallest peak-to-trough 
wave height of interest. For the present SEASAT application this corresponds 
to an antenna smaller than 1.1 meter and a pulse width of about 3 nsec, which 
is readily achievable in the recommended design. 

Note that in the comparison of the two modes of operation the 
emphasis is placed on the ability to faithfully reproduce the wave height 
density or distribution function so that an estimate of wave height may be 
obtained. It is also seen that both modes required a wide bandwidth signal 
(—330 MHz) in order to resolve the smallest (1 meter) wave height of interest. 

As will be shown, this bandwidth requirement is basic to the measurement of 
wave height regardless of the technique employed. 

Consider now the two frequency correlation techniques suggested by 
Weissman^ for the measurement of wave height. With this method two very long 
(essentially CW) signals separated in frequency by the amount Af = - f^ 

are simultaneously transmitted. Two homodyne receivers are then used to 
separate the return signals at frequencies and and measure the cross 
correlation R(Af) at the frequency separation Af. Thus changing the transmitter 
frequency separation Af provides a measure of R(Af ) . It can be shown ^ that 
the cross correlation, R(Af), is simply the Fourier transform of the radar 
cross section versus range CT(R). Now if the radar is operating in a narrow 
beam mode (for which the antenna requirement for has already been 

shown to be impractical at satellite altitudes) the Fourier transform of 
0 ( R) , FT[a(R)], is the characteristic function of the specular point density. 

Note also that the wave height resolution of this implementation is the reciprocal 
of the transmitter frequency separation, 1/Af, and hence would require the same 
330 MHz for the measurement of 1 meter (peak- to-trough) wave heights as the 
short pulse radar. The point to be made here is that the measurement of the 
wave height density or characteristic function both require the same overall 
system bandwidth and antenna requirements. 
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Finally, consider operating the two frequency interferometer in a 
wide beam mode for measuring wave height. Here, the measured output R(Af) 
would again correspond to the FT[cr(R)], but in this situation the transform 
is not only that of the cumulative wave height distribution (corresponding to 
the leading edge of a(R)) but also includes the so called plateau region of 
a(R). Unfortunately, most of the energy in a(R) is contained within the 
plateau region and hence also in FtCct(R)]. Thus, changes in wave height 
would have little effect on the main response of FT[<j(R)] and would probably 
have to be determined by the measurement of the sidelobe structure of FT[o(R)]. 
While such a measurement technique has not been examined in depth (and in 
fact might even be feasible) it certainly cannot compare favorably with the 
short pulse wide beam radar which provides altimetry as well as wave height 
information. Table 2.1 summarizes the more important characteristics of a 
short pulse radar and two frequency interferometer technique for measuring 
wave height. 

From the above discussion it is seen that a short pulse wide beam 
radar is the most appropriate (and recommended) technique for measuring both 
altitude and wave height. The following sections summarize the design of such 
a radar, which includes an on board MLE processor. 

2 . 2 System Design Summary 

The design of a high resolution satellite radar for measuring altitude 
and wave height is described on the following pages. The resulting design achieves 
all specified performance requirements. These performance requirements are 
summarized in Table 2.2 and the system design is summarized in Table 2.3. 

A simplified block diagram of the proposed radar is shown in 
Figure 2.2. The CHIRP transmitter generates a 2.8 psec 360 MHz bandwidth 
linear FM signal. When the chirp transmitter receives both a PRF trigger 
and a transmit pulse gate, the output signal is fed to a TWT string identical, 
to that used in GEOS-C. Full deramp pulse compression is achieved b;y 
generating a second identical linear FM signal by means of a properly timed 
range trigger from the MLE processor. This second linear FM signal is not 


1 - 2-6 



Table 2.1 


Short Summary of Wave Height Measurements 


Quantity 

Type of Radar 

Short 

Pulse 

Two-Freq. 

Interferometer 

Resolution 

1/B 

l/Af 

Max date rate 

B 

1/ T 

Energy on 
target 

T 

t( 50% loss) 

Measured 

output 

a(R) 

FT [a (R) ] 

a) Narrow 
beam 

Density 

Characteristic 

b) Wide 
beam 

Cumulative 
dis t . 
function 

F.T.Cc.D.F. + 
Plateau] 


Table 2.2 

System Performance Requirements 


I. 

Geodetic accuracy 

50 cm 

II. 

Topographic resolution 

10 cm RMS (7 cm allocated 



to system error) 

Ill, 

Wave height range: 

1-20 m crest-to-trough 


Accuracy: 

max (0.5 m.,25%) 

IV. 

Correlation between 
pulses 

' <l/e 

V. 

Oceanographic phenomena 
of interest (maximum 
spacial frequency) 

0.25 Hz 
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Table 2.3 


Design Summary 


I . 

Orbit Parameters . 



a) Height 

725 km 


b) Inclination 

82° retrograde. 


c) Eccentricity 

0.0064 maximum 

II. 

Radar Parameters 



a) Antenna beamwidth 

2.6° 


b) Pointing accuracy 

20 q =0.1° 


c) Antenna gain 

36.4 dB 


d) Peak power 

2 kW 


e) System Losses other than 
processing losses in pulse 
compressor 

5 dB ' 


f) Noise figure 

5.5 dB 


g) Frequency 

13.9 GHz 


h) Uncompressed pulse width 

2.8 M-s 


i) Uncompressed pulse bandwidth 

360 MHz 


j) Compressed pulse width 

3.0 ns 


k) Compression ratio 

1000/1 


1) PRF (uncorrelated returns) 

max x 

1.6 kHz 


m) PRF 

1500 Hz 


n) S/N (single pulse) 

10 dB 


o) Ocean cross section 

+6 dB 


p) Receiver weighting* 

-26 dB Modified Taylor 

1 

q) Pulse compression * 
processing loss 

0.55 dB 


r) Main lobe broadening due to 
tapering* 

237. 

III. 

Pulse Compression 



Type 

Full deramp stretch , 


Range processing 

Analog filter bank 


Filter bank 

Discrete passive 

1 

Number of filters 

60 

1 

} 

Frequency range 

9.2 to 20.8 MHz 

( 

Filter bandwidth 

385 kHz | 

1 

Output data form 

Two TTL parallel words 

A, Range bin number, 6 bits 

B. Range bin amplitude, j 

13 bits 

Included in design but considered optional; 

| 

! 



Table 2.3 (cont.) 


III. 

Pulse Compression 



Time required for full 

450 microseconds, max. 


sampling 



A/D sampling frequency 

1 MHz 

IV. 

Linear FM Generation 



Type* 

Surface wave 


Bandwidth 

60 MHz 


Multiplier chain 

X6 


Pulse length 

2.8 M>s 


Linearity of FM 

0.2% 


Peak frequency deviation 
(one circle of variation 
across pulse) 

25 kHz 

V. 

Altitude and Wave Height Processing 



Type 

MLE processor 


Implementation 

Digital interface with on- 
board computer 


Interface data rate 

200-360 Hz 


Output data rate 

3-5 words @ 1-10 per second 


Tracking bandwidth 

1.0 Hz 
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Figure 2.2 High Resolution Radar Altimeter Block Diagram 



transmitted but instead fed to the deramp receiver where it is mixed with the 
ocean return signal* This process converts the return signal as a function 
of time into an identical signal as a function of frequency. 


This return signal as a function of frequency is then fed at an IF 
frequency into the range processor unit. There an analog filter bank con- 
sisting of 60 discrete filters is used to provide 60 samples of the leading 
edge of the return signal (i.e., cumulative wave height distribution function). 
The 60 range samples are then detected and fed into an analog multiplexer and 
A/D converter. The output of the range processor unit then consists of two 
TTL parallel words, range bin number (6 bits) and range bin amplitude (13 bits). 


A more detailed breakdown of the recommended design for the chirp 
transmitter, deramp receiver, and range processor is shown in Figures 2.3, 

2.4 and 2.5 respectively. This particular design is the result of a previous 
NASA/Wallops study An in depth description of the rational and tradeoffs 
involved in arriving at this design is available in the above referenced final 
report. The most important aspect of the above design is the selection of a 
full deramp pulse compression technique followed by an analog filter bank to 
separate individual range returns. With this technique the A/D converter 
bandwidth is less than 1 MHz as compared to other forms of pulse compression 
which would require a rather impractical 330 MHz A/D converter. 


The MLE processor accepts the digital output of the range processor 
and simultaneously performs an optimum (minimum variance) estimate of epoch, 
wave height, and signal- to-noise ratio. In addition, the MLE processor also 
provides a range trigger to the chirp transmitter unit for deramping the return 
signal and hence closing the altitude tracking loop. As shown in Figure 2.6, 
the digital input from the range processing unit is fed into a buffer storage 
and inverted to produce negative digital video samples, -V, which when summed 

A A 


with the estimated mean power return V produces a difference signal, V-V. 

This difference signal is then normalized by dividing by the variance of the 
—2 

estimated return V , multiplied by the appropriate parameter weighting function 


9? o 


, or ) , and, range summed in each loop to produce the error signals 

aa 0 30 h - ! 
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Figure 2*6 Onboard MLE Processor 









el , , and el * After passing through the loop filters, the estimates 

”o o °h 

are fed back as follows: the epoch estimate is used to control the range 

trigger of the chirp transmitter for deramping the ocean return signal. 

The wave height estimate is used to control the table look up of the mean 
return signal as well as the parameter weighting functions* The signal- 
to-noise ratio estimate is used as a scale factor in the generation of the 
estimated mean return and as a scale factor in the generation of the epoch 
and wave height weighting functions. 


The implementation of a MLE processor could take one of two forms; 
a complete special purpose digital logic unit or a digital interface with a 
minicomputer. If a minicomputer is included in the proposed SEASAT-A con- 
figuration, a digital interface would be the recommended approach* The 
interface data rate would be 62 words at a nominal 200 - 360 Hz rate* The 
MLE output data rate would be 3 - 5 words at a 1 - 10 per second rate* 


Complete details of the mathematical development of the MLP pro- 
cessor are given in Reference (1), and further delineated in Section 3,0 
and Appendix A of this report* Again, it should be pointed out that some 
sort of adaptive processor — whether it be an MLE, MMSE, or ASG — is 
necessary (and recommended) to achieve the performance requirements on board 
the spacecraft. A comparison of tracking accuracies, Figure 2.7, shows that 
a 5*0 meter RMS wave height, a factor of four reduction in RMS tracking error 
is achieved by the MLE when compared with a conventional half power split- 
gate tracker. The accuracies and biases associated with these split-gate 
trackers (conventional and adaptive) are examined in detail in Section 4,0 
of this report. 
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Figure 2.7 Comparison of Tracking Accuracies of Split-Gate Trackers and Adaptive Processors 
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3.0 


MAXIMUM LIKELIHOOD ESTIMATOR 


A joint Maximum Likelihood Estimator (MLE) which would simulta- 
neously estimate the epoch, wave height, and signal-to-noise ratio for a 
satellite altimeter was originally described in Reference (1) . This sec- 
tion of the report describes a method of implementing a MLE processor on a 
digital computer. Estimation accuracies as well as biases associated with 
range sidelobes are given for typical SEASAT-A altimeter parameters. The 
intention of this section of the report is to describe results and explain 
the simplicity of the concepts involved in the implementation. The theo- 
retical basis for all of the material described here is provided in Appendix 
A. 

Section 3.1 describes a feedback loop implementation of the MLE 
processor. It is shown that for such an implementation, the logarithms of 
the likelihood ratio must converge with time to a local minimum of the 
likelihood function. A step-by-step description of the implementation is 
presented. Simulated examples which illustrate the shape of the functions 
used to generate the error signals of the processor are described in detail. 
Some of these functions are shown to be generalizations of the gating func- 
tions used in split-gate trackers, thus providing further insight into the 
nature of the feedback loop implementation. 

Section 3.2 describes the approximate "size 11 of the MLE processor 
when implemented as a software program in a general purpose digital computer. 
It is shown that the total data storage required is about 2000 words, and 
the total program storage is also about 2000 words. These results indicate 
that a relatively modest mini-computer would be adequate for implementing 
the processor. 

Section 3.3 presents the theoretical estimation accuracies of the 
MLE for typical SEASAT-A system parameters. An evaluation example is given 
for each of the three parameter estimates. A table of estimation accura- 
cies shows that the performance requirements for SEASAT-A could be achieved 
with a signal-to-noise ratio greater than 5 dB even at 20 meter significant 
wave height. 
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Section 3.4 provides expressions for the biases in the three param- 
eter estimates caused by range sidelobes of the pulse compression network. It 
is shown that all biases can be maintained at a negligible level (less than 
1 percent) , provided that the average RMS sidelobes are less than or equal to 
-50 dB. Since these sidelobe levels can be achieved in practice, the result- 
ing biases for the MLE should not cause any problems. 


3. 1 Feedback Loop Implementation 

The implementation of a joint Maximum Likelihood Estimator (MLE) 

for epoch, normalized surface reflectivity, and RMS wave height (T , a°, O ) 

o h 

is examined in this section of the report* Here, the intention is to explain 
the simplicity of the concepts involved, at the expense of rigorous theoret- 
ical justification. The theoretical developments are given in Appendix A 
of this report. 


First consider the recommended altimeter design, which consists 
of a full deramp pulse compression technique followed by an analog filter 
bank to separate individual range returns. Then, the altimeter receiver 
can be modeled as (nearly) matched filtering, square law detection and 
sampling as shown in Figure 3.1. The description of the MLE and its pro- 
perties are simplified with this model, since the problem is now constrained 
to that of finding the optimum processor for the sampled video outputs. In 
the sense that only the video processing is optimized, the resulting pro- 
cessor is subop timal ; however, this is not thought to be a serious limita- 
tion. Aside from the fact that the optimum video processing to be described 
will show considerable improvement over existing tracking schemes, it also 
readily lends itself to a practical implementation with a small scale digital 
computer. 


The procedure to be followed will be to find a suitable approxima- 
tion to the likelihood function (i.e. , the joint probability density) of the 

video outputs given the three parameters of interest T , o° and a . Then, 

° h 

for any set of observed video samples, the joint MLE are those values (*f , 
o o 

a , a h ) which maximize the logarithm of the likelihood function. The max- 
imum, of course, occurs at a point in the parameter space for which the 
three partial derivatives with respect to the parameters are zero. This 
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Figure 3.1 

Functional Model of the Altimeter Receiver 
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fact suggests a technique for implementing the MLE which uses the partials 
as input error signals to three feedback loops which, in turn, solve the 
maximizing conditions by forcing the error signals to zero. 

The form of the likelihood function for the sampled video outputs 
of Figure 3.1 is simple, if it is assumed that the samples are separated by 
at least one range resolution cell. In this case, the correlation between 
samples is negligible, and since the underlying process is Gaussian, it is 
reasonable to assume them to be independent. Thus, the likelihood function 
has the form: 


A (V 


V 

11 ’ 12 


V . /t , 
nk o 


o 

O , 


V - " V 1 

i,k 


exp( " V ik /V i } 


where: 


V ik 


th 

is the sampled video output from the i — 
range cell on the kit pulse. 


and 


V. 

i 


is the expected value of the video from 
the iit range cell. 


(3.1) 


That is, the square law detected outputs form an independent 
exponential process. The dependence of the likelihood function on the 
three parameters t , a° , is contained entirely in the variation with 

range of the mean value V^, Thus, taking the logarithm of (3.1), and 
differentiating, the MLE estimates must satisfy the three equations. 


0 = 


3 log A 

3 a 


-2 


= - £ 
k 


(V i 


V ik> 


3 V. 

1 

S a 


(3.2) 


where "a" stands for any of the three parameters T , a° and o^. 

The problem of determining the joint MLE of the parameters requires 
solving Equations (3.2). One technique for achieving this is to use the 
negative of the partial derivatives in Equations (3* 2) as inputs to integrating- 
feedback filters to derive the MLE estimates as shown in Figure 3.2* 
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That is, since the estimates are derived from integrating filters, 


one has: 


dT 
o 

dt 


K5 log A 

a t 


d q _ _ K d log A 


d t 


aO 

9 a 


(3.3) 


d J h 
d t 


_ K 9 log A 


a 


where K is a gain factor associated with the integrators. 

Now, the time rate of change of the logarithm of the likelihood 
function is: 

^ Ae* + + £ < 3 - 4 > 

o do h 

or, from Equations (3.3), 


d log A 
dt 


- K {( ^ °a S - ) 2 + ( ~ 1? &- ) 2 + ( dg^M ) 2 ] 
o 9 a h 

0 


(3.5) 


Thus, as a function of time, the logarithm of the likelihood 
ratio must decrease and converge to a local minimum of the likelihood func- 
tion. If the initial parameter estimates are close to the correct ones, 
then the loops in Figure 3.2 will converge to the joint MLE of T , q°, and 

V 

The practicality of implementing the MLE depends on the complexity 
of the partial derivative calculator in Figure 3.2. The steps required to 
compute the three partial derivatives are diagrammed in Figure 3.3. As can 
be seen from that figure, if the mean video versus range plus the three 
partial derivatives are available, then the computations required to derive 
the error signals are relatively minor. 
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Figure 3.3 Generation of the Error Signals for the MLE 
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Figure 3,4 Onboard MLE Processor 



Figures 3.2 and 3,3 give a fairly pleasing intuitive picture of 
the operation of the MLE processor. In Figure 3.3, it is seen that the 
error signals for the tracking loops in Figure 3.2 are derived as follows. 
First, an error signal versus range, V - V, is generated which represents 
the difference between the measured return as a function of range and the 
estimated mean return. Secondly, this error signal is weighted inversely 
proportional to its variance. That is: 


E (¥- V) 2 = V 2 (3.6) 

Thus, after normalization, the signal may be loosely described as having 
uniform information content. Finally, the normalized error-versus-range 
signal is "gated" by the partial derivative of V, with respect to the 
parameter of interest, and summed over range. The effect of multiplying 
(gating) by the partial derivative is to emphasize those range bins which 
are most affected by variations in the parameter of interest. 


Thus, a very simple step-by-step description of a feedback loop 
implementation, Figure 3.4, for solving the minimizing conditions 


0 



£ 

i 


< v - V £ 

^2 da 


and achieving the joint MLE is as follows: 


Step 1. Compare an estimate of the mean return signal with the sampled 

data to obtain a difference signal V - V,, * 

rk 

Step 2. Normalize the fluctuation of this difference signal by dividing 
it by the estimated variance of the return V . 


Step 3. Weight this normalized difference signal to maximize its sensi- 
tivity to the parameter being estimated (multiply by av/da). 

Step 4. Integrate the resulting signal over range to obtain the error 
signals e . 

Step 5. Integrate the error signal e_ with time to obtain the parameter 
estimates a. 

Step 6. Use the parameter estimates & to obtain a new estimate of the mean 
return signal V so as to drive the error signal to zero. 
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3*1.1 


Simulated Examples 


The nature of the error signals generated by the feedback loop 
implementation of the MLE processor becomes more apparent from an examina- 
tion of several simulated examples. These examples illustrate the shape 
of the functions utilized in generating the error signals. Three examples 
are given, one for an error in epoch only, one for an error in wave height 
only, and one for an error in ocean reflectivity only. The parameters 
used in simulating these examples are as follows: 

Data: a° = 10 dB, ^ = 2m, T 0 = 0 (true values) 

Sample Interval = 1,0 m 
Number of Samples Averaged = 5 

Case I: ¥ error = 1.0 m 

o 

No error in a or 

h 

Case II: 6, error = 1.0 m 

n 

No error in a or t 

o 

Case III: 6° error « 3.0 dB 

No error in a , or t 
h o 

The shape of the functions utilized in generating the error sig- 
nals for Case I,epoch error only, are shown in Figure 3.5-a through 3.5-f. 
Figure 3.5-a shows the estimated mean V, the data samples V^, and the true 
mean of the data samples V. Here, and throughout all the examples to be 
described, each data sample V. is the average of five independent samples. 
That is, it is assumed that some pre-processing in the form of pulse-to- 
pulse averaging is performed prior to inputting the data to the MLE pro- 
cessor. In the actual MLE implementation, the amount of pre-processing 
(averaging) should be such that the averaged sample values are input to 
the processor at a maximum rate of 10 per second. 

In addition to showing the 1.0 meter displacement between the 
true and estimated epoch, Figure 3,5-a also clearly illustrates the fact 
that the variance of the data samples is proportional to the mean value as 
shown by the fact that the fluctuation of the data samples increases as the 
mean value increases. This effect is even more pronounced in Figure 3.5-b, 
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Figure 3.5 Epoch Error Example (Case I) 
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which shows the difference signal between the estimated mean and the data 

a - ■ 

samples, V-V., as well as the difference between the estimated mean and 
a i 

true mean, V-V. When this difference signal is normalized by multiplying 

— 2 

by the reciprocal of the estimated variance 1/V shown in Figure 3.5-c, 
the normalized difference signal shown in Figure 3,5-d is obtained. The 
mean value of this normalized signal is also shown in the figure. Note that 
the fluctuation of the data samples of the normalized difference signal has 
been reduced to about +0.5, as compared with the difference signal shown 
in Figure 3.5-b, which has fluctuations on the order of + 5.0. This type 
of normalization of the data for reducing and equalizing the fluctuation 
is characteristic of all maximum likelihood estimators. 

When the normalized difference signal is multiplied by the epoch 
weight shown in Figure 3.5-e, the signal shown in Figure 3.5-f is obtained. 
Here, several observations can be made. The shape of the epoch weight shows 
that the maximum weight is placed at the mean of the estimated return sig- 
nal, Furthermore, since this weight is the partial derivative of the esti- 
mated return, it is in fact the estimated wave height probability density 
function. Thus, the width of this weighting function is equal to the esti- 
mated RMS wave height cr^* which of course, changes width with sea state. 

This epoch weight of the MLE processor is analogous to the early gate of a 
split-gate tracker. It is, of course, a more general type early gate in 
the sense that the weighting is continuous, and the width changes with wave 
height. 

If the epoch weighted normalized difference signal shown in Figure 
3.5-f were integrated (summed) over range, a single value of epoch error 
signal would be obtained. Integrating successive samples of in the 

epoch loop filter would thus provide a new epoch estimate T . While it 
would have been desirable, time and funding limitations prevented includ- 
ing the shape of the error signals as a function of error for this as well 
as the other simulated examples. It may be shown, however, that all the 
error signals are strongly coupled. That is, an error in epoch estimation 
only (as in this example) produces not only an epoch error signal, but also 
error signals in the reflectivity and wave height feedback loops. The same 
sort of thing happens when there is an error in estimating either wave height 
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or reflectivity. Whether or not this coupling of the error signals pre- 
sents a problem for the feedback loop implementation of the MLE remains 
to be determined. Since the coupled error signals are such that the total 
error in all three loops is driven to zero as rapidly as possible, there is 
reason to believe that the coupling will not be a problem. On the other 
hand, if error coupling proves to be a problem, the situation can be cor- 
rected by means of a decoupling network. Such a network would entail esti- 
mating (and inverting) the error correlation matrix R(0), which would add 
to the complexity of the processor. To resolve these questions concerning 
the effects of error coupling requires a detailed study of the eigenvalues 
and eigenvectors of the error correlation matrix which was beyond the scope 
of the present program. 

A description of the shape of the functions utilized in generat- 
ing the error signals for Case II, wave height error only, shown in Figure 

3.6- a through 3.6-f, and Case III, reflectivity error only, shown in Figure 

3.7- a through 3.7-f, is essentially the same as the previous example except 
for the weighting functions. An examination of the wave height weight, 
Figure 3.6-e, shows a maximum negative weight applied one sigma (RMS wave 
height a^) prior to the estimated epoch and a maximum positive weight 
applied one sigma after the estimated epoch. Functionally, this weight is 
formed by the product of a linear ramp and the estimated wave height density 
function. A gating arrangement analogous to this wave height weight would 
entail centering two gates symmetrically about the half power point of the 
return signal. The difference between these two gates would thus provide 

an estimate of the mean slope of the leading edge of the returns signal 
and hence a measure of wave height. (This is, in fact, the method sugges- 
ted by MacArthur , Reference (2), for estimating wave height.) The wave 
height weight of the MLE processor is of course, more general in the sense 
that the weighting is continuous, and the separation and width of the posi- 
tive and negative portions change continuously with wave height. 

An examination of the reflectivity weight, Figure 3.7-e, shows 
the maximum weight being applied well after the estimated epoch. Since 
this weight is the partial derivative of the estimated return with respect 
to a , it is simply a scaled version of. the estimated return, and hence, 
corresponds to the estimated cumulative wave height distribution function. 
Quite obviously, this reflectivity weight is the more general analogy to 
the late gate of a split-gate tracker. 
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c) Normalization Based on Estimated Variance 


Figure 3*6 Wave Height Error Example (Case II) 
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In summary then, it is seen that the feedback loop implementation 
of the MLE processor is simply a generalization of the split-gate tracker 
and gating technique for estimating wave height* (Actually, the gating 
technique for estimating wave height was developed as a simplification of 
the MLE wave height weighting function.) Aside from the more general weight- 
ing functions, the feedback loop implementation of the MLE also differs from 
the gating techniques, in that a replica (or template) is compared with the 
return signal, and the fluctuation of the data is then reduced by normaliz- 
ing with the estimated variance. It might also be noted that if the nor- 
malization by the estimated variance is eliminated, the feedback loop 
implementation would be that of a MMSE processor! 

3.2 Computer Sizing 

In this section, an approximate "size" of an MLE is obtained. 

For this estimate, it is assumed that the processor will be implemented 
as a software program in a general purpose computer. The size is then 
determined as the total program and data storage required, plus the total 
compute time required during each data cycle. The results indicate that 
a relatively modest mini-computer is required. It is shown that the total 
data storage is approximately 2000 words and the total program storage is 
about 2000 words; thus, only about 4 k memory words are required. Further, 
assuming a compute speed roughly equivalent to the NOVA, the total comput- 
ing time per data cycle is about 28ms. If the data is averaged for 0.1 
sec (100ms) before the MLE acquires it, then the computer will be rela- 
tively lightly loaded, 

3*2.1 Approach 

This study was intended only as a 11 quick look 11 estimate; thus, 
the approach was correspondingly simple. First, a detailed block diagram 
of the processor was developed. This is shown as Figure 3.8. In that 
figure, the data is supplied to the MLE at 10 Hz. Thus for each .1 sec 
data cycle 62 power values are transferred into the buffer. The levels 
represent the integrated power for .1 sec (150 pulses) for each of the 
60 range gates on the leading edge plus a noise gate and a plateau gate. 

The diagram to the right of the dotted line represents a detailed break- 
down of the computations performed in the MLE. 


1-3-19 



Ui 

N> 

O 



Figure 3*8 MLE Implementation Block Diagram 








The computation speed estimate was obtained from Figure 3.8 as 

follows: 

• On a per "data cycle" basis, count the number of arithmetic 
operations and I/O requests indicated in the diagram. 

• Using "typical" computer speeds, and the above counts, 

. estimate the compute time per data cycle. 

• Apply 25 percent increase for program housekeeping, (loop 
initialization, data shuffling, etc.) 

• Apply 15 percent increase for the real time executive. 

The above approach is obviously somewhat crude, however, to obtain 
better estimates, one would have to write a sample code and then count the 
number of operations in the code. Further, that approach would probably not 
yield particularly accurate estimates until it was actually run in a real- 
time environment. Either of these approaches are beyond the scope of this 
study. 

The percentage increase in the running time for housekeeping and 
the real-time executive have not been documented* They are "engineering 
judgments" based on discussions with several people who have had experience 
in programming mini-computers in a real-time environment. 

The computation of the core storage requirement was obtained from 
Figure 3.8 as follows: 

• Data Storage 

From the block diagram, count the data storage. 

• Program Storage 

For the tracking phase, identify major function blocks 

and estimate the number of instructions per block. 

Assume that the code for the acquisition phase approx- 
imately equals that for the tracking phase. 

Add 250 instructions for initialization. 

Add 200 instructions for the executive. 

Again, the constants for the initialization and the executive 
represent "engineering judgment" rather than documented fact. 
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3.2.2 


Details of the Estimates 


As stated previously, it is assumed that 60 leading edge samples 
plus one noise gate plus one plateau gate are transferred to the MLE each 
data cycle. The data cycle rate is 10 Hz; thus, for a 1500 Hz PRF, the 
altimeter averages 150 pulses before transferring the data to the MLE. 

The total number of arithmetic operations per data cycle is given 
in Table 3.1. To convert the totals to required computing time, it is 
assumed that, on the average, about five total instructions requiring approx- 
imately 10 p,sec are needed for each addition or subtraction. For multiplica- 
tion and division, about eight instructions and 20 ^sec are required. These 
numbers assume that the registers are loaded from memory, and the results may 
be normalized before storage. The compute time corresponse approximately to 
fixed point operations on the NOVA II computer built by Data General Corporation. 

Table 3.2 summarizes the total compute time per data cycle. The 
compute time, housekeeping, and executive have been explained previously. 

The I/O handling represents a pessimistic estimate. It is assumed that the 
computer supplies trigger information to the altimeter at the PRF, and thus 
must provide 150 numbers per data cycle. Further, it is assumed that for 
each number, the system must process an interrupt and store working regis- 
ters before processing the data. This is assumed to take 30 p.sec per 
interrupt. 


From Table 3.2, it is seen that the total compute time is esti- 
mated to be ^ 28 ms. Thus, the computer would be 28 percent loaded. This 
seems to be an adequate margin for a computer dedicated entirely to the 
MLE function. Conversely, this probably represents a significant load on 
a time share system in which the computer would be performing additional 
system functions. 

The data storage and program storage are detailed in Tables 3.3 
and 3.4. The numbers there should be self-explanatory. The major data 
storage is in the tables of the wave height cumulative distribution func- 
tion (CDF) and the probability density function (PDF). These tables are 
used to generate an estimate of the mean return, V. 
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TABLE 3-1 


Speed Requirements 


FUNCTION 

ADDS 

MULTIPLIES 

DIVIDES 

Noise Integration 
and Normalization 

1 

1 


Plateau Normalization 


1 


Error Signal Generation 




V 

62 

62 


e 

o 

62 

62 

62 


62 

62 


e a 

o 

62 

62 


\ 

62 

124 


Loop Filters 




• 

T o 

2 

3 


T 

0 

150 



a 

o 

3 

3 


a h 

3 

3 


Calib. Normalization 


4 

1 

Ramp Generator 

120 

2 

2 

Offset 

150 



Totals 

739 

389 

65 
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TABLE 3.2 


Compute Time Per Data Cycle 


"Typical" Computer: 




Add Time 

5 Instructions 

10 

\isec 

Multi/Div Time 

8 Instructions 

20 

\isec 

Compute Time 
Per Data Cycle 

16.5 msec 



Housekeeping ( 25%) 

4.1 msec 



Exec. (157.) 

2.5 msec 



I/O Data Handling 

4.5 msec 



Total 

27.6 msec 
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TABLE 3.3 


Data Storage Locations 


Input Buffer 

124 


T Buffer 
0 

150 


Output Data 

6 


Intermediate 

20 


Storage 


Loop Filters 

9 

(Constants) 


6 

(Previous Values) 

Ramp Generator 

6 


Calibration 

20 


CDF Table 

600 

at 5 cm Resolution 

PDF Table 

600 

for 20 m Wave Height 


< 2 k 
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TABLE 3.4 


Program Storage 


Tracking Mode: 


Major Block 

Number of Instructions 

Input Buffer 

50 

Noise Integrator 

75 

Error Signals 

200 

Loop Filters 

225 

Ramp Generator 

50 

Calibration 

25 

Total 

625 

Acquisition Mode: 

625 

Initialization 

250 

Executive 

200 


2 k 
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The program storage estimates again were in the nature of Engi- 
neering judgments” rather than actual instruction counts* It was beyond 
the scope of this task to do the coding which would be required for more 
accurate counts., It should also be pointed out that these counts assume 
that all code has been generated from scratch. In particular, no manufac- 
turer supplied standard I/O, or executive routines have been used. Such 
routines would usually be more general purpose than is needed, and would 
increase the storage requirements considerably (e.g., as much as a factor 
of 2 or more) . 

In sum, the size requirements from the first cut estimates appear 
to be modest, but not trivial. Approximately 4k of memory is required, and 
a computing speed comparable to existing mini-computers (NOVA II) appears to 
be adequate with a reasonable safety margin (28 percent loaded). This is 
true, provided the computer is dedicated to the MLE, or performs only a few 
additional tasks. 

3*3 Estimation Accuracies 

The theoretical accuracies of the joint maximum likelihood esti- 
mates of epoch, wave height, and reflectivity are evaluated in this section 
of the report. These accuracies are given in terms of the standard devi- 
ation of the estimates due to random fluctuations only. System errors 
(such as clock and timing errors) as well as errors associated with the 
feedback loop implementation as described in Appendix A, are not included 
in the evaluation. The theoretical basis for the expressions used in eval- 
uating these accuracies was originally derived in Reference (1), and is 
rederived and expressed in a much simpler and more useful format in Appen- 
dix A. An example showing the method of evaluation is given for each of 
the three estimates. A table of estimation accuracies as a function of 
signal- to-noise ratio and wave height is provided for each parameter esti- 
mate under typical SEASAT-A conditions. 

The following typical SEASAT-A parameter values are used in each 
of the sample calculations: 
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Number of Pulses 


1500 


Range Resolution 

Data Interval 
RMS Wave height 
Signal-to-Noise Ratio 


N 

cAt 

2 


cT 


max 


a 


.5 m 

23 m 
5 m 
10 dB 


Here the data interval is assumed to contain half of the 60 samples on the 
leading edge of the return plus an additional 16 samples contained in a 
late gate. This results in a total of 46 independent range samples per 
pulse, 

3,3,1 Epoch Example 

From Appendix A the variance of the epoch estimate is 

a T T 2 = [if 1 R(e) _1 ] T T (3.7) 

o o o o 


where 


of c " 1 is 


R(e )" 1 = apA t d" l c -1 d -1 

.-1 


(A-67) 


= diag(a, (3 , (3) is a diagonal matrix, the (i,j) element 


and 



l + dF n , li lj 


d ' - « <TT1) 2 


(A-66) 


(A-68) 


Table A- 1 of Appendix A tabulates F„ versus signal-to-noise ratio from 
-10 dB to 30 dB in 5 dB steps. The parameter f3 is related to RMS wave 
height via 



(3.8) 


1-3-28 



and the constant a= .3227 arises from fitting a linear ramp to the leading 
edge of the return signal in a minimum mean square sense. Also, the sub- 
scripts of Equation (A-66) are related to the parameter being estimated via: 


1- **»a 

2- ►t 

c 

3- ~-P 

Now substituting in Equation A-68, 


(3.9) 


, , a T max . a .2 ,.3227 , w 10v2 

j = (——5 — ixrrp ■ (23) - 

= .4003 

From Table A-I (using 10 dB signal-to-noise ratio), 

= 2.42, F^2 = 3.75, F 22 = 9,82, and from Equation A-66 


T T 

c ° ° = C 22 = F 


d F 2 


22 1 + d F u 12 


= 9.82 - 


.4003 


1 + (.4003) (2.42) 


(3.75)' 


= 6.961 

Substituting Equations A-67 into Equation 3,7 yields: 


T T 

o* T = (a^At) N -1 ff 2 C ° ° 
o o 


or 


°R 




fair—) a h n ’ 1 C T ° T ° 


(3.10) 


• 3227 (.5) (5)0^) 6.961 


= .061 meters or 


6.1 cm 
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3,3.2 Wave Height Example 

The variance in estimating the parameter (3 

a p/ = [N_1 R(e) " 1] pp = (ct P At) rl P 2 C ^ 

after substituting Equation A-67. 


(3.11) 


From Table A-I, F 1n = 2.42; = -2.42; F,^ = 5.14; and as in the 


11 

previous example d = .4003, Thus: 


- C 33 - F 


d F. 2 


33 


33 1 + d F 1X 13 


= 5.14 


.4003 


1 + (.4003) (2.42) 


(2.42) 


= 3.95 

Converting from the parameter (3 to and using Equation 3.11 


2 2 / , cAt s M - 1 

' \J a < -— ) “h" 

- V ,3227 i^ (5> w 


.pp 


= *046 meters or 


4*6 cm 


( 3 . 12 ) 


(3*13) 


3.3.3 Reflectivity Example 

The accuracy with which ocean reflectivity can be estimated is 
given in terns of the variance of the signal-to-noise ratio 

O aa 2 = [n" 1 R(e) -1 ] aa = (00 At) N -1 a 2 C aa (3.14) 

where again the substitution of Equation A-67 has been utilized for the 
inverse of the error correlation matrix* 
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From Table A-I, F n = 2.42; and d - .4003. Then: 


„aa „11 
C = C = F 


d F. 2 


11 1 + d F X1 11 


- 2.42 


* 4003 (2.42) 2 


1 + (.4003)(2.42) 


- 1.229 


and 


°aa = a \/ a ^ N_1 <?* 


= 10 


,L 3227. 0 . 5) (. 2X [ — 

V 1500 * y 


.051 


(3.15) 


Estimation accuracies versus signal-to-noise ratio and significant 
wave height are shown in Table 3.5 for typical SEASAT-A system parameters. An 
examination of the table entries shows that even at 20 meter significant wave 

height, the performance requirements for SEASAT-A could be achieved with a 
joint MLE operating at signal-to-noise ratios greater than 5 dB. This assumes 
that for the altitude estimate, the system errors (clock, timing, etc.) can be 
maintained to within 7.0 cm, and that any errors associated with the feedback 
loop implementation of the MLE processor do not significantly change the theo- 
retical performance values given in Table 3.5. 

3.4 Range Sidelobe Bias 

The biases associated with range sidelobes (due to pulse compres- 
sion) that have not been included in the model for the estimated mean power 
return are summarized in this section of the report. The rather complex 
theoretical developments and evaluation of these biases are provided in 
Appendix A. There, it is shown that biases in estimating altitude, wave 
height, and signal-to-noise ratio caused by range sidelobes can be expressed 
in the following simple form: 
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TABLE 3.5 

Estimation Accuracies of a Joint MLE 


SYSTEM PARAMETERS 


PRF 

= 1500 Hz 

Loop Bandwidth 

= 1 Hz 

Resolution 

= .5 m 

# Range Samples 

= 60 

# Samples in Late Gate 

= 16 


STANDARD 

DEVIATION IN ALTITUDE (cm) 


Significant Wave Height (m) 

S/N (dB) 

5 

10 

20 

0 

5.7 

8.4 

13.1 

5 

3.2 

4.8 

7.8 

10 

2.5 

3.7 

6.1 

20 

2.1 

3.1 

5.3 

STANDARD DEVIATION 

IN WAVE 

HEIGHT (cm) 


Significant Wave Height (m) 

S/N (dB) 

5 

10 

20 

0 

5.9 

8.5 

12.6 

5 

3.0 

4.4 

6.7 

10 

2.0 

2.9 

4.6 

20 

1.4 

2.1 

3.5 

STANDARD DEVIATION IN 

SIGNAL- 

TO-NOISE RATIO 

(ABSOLUTE UNITS fNot dBl ) 


Significant Wave Height (m) 

S/N (dB) 

5 

10 

20 

0 

.008 

.008 

.009 

5 

.017 

.017 

.019 

10 

.044 

.046 

.051 

20 

.402 

.422 

.472 
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AR = 880 

"h * 600 °h n 

Aa = 600 a y^ (3.16) 

where y^ is the average RMS sidelobe level. In deriving these expressions 
it has been assumed that the 3 dB (one way) antenna beamwidth of the altim- 
eter is 5° and the satellite altitude is 725 km. A rather surprising fea- 
ture of these results is that the biases are independent of signal-to-noise 
ratio! Furthermore, the biases are all proportional to average RMS side- 
lobe level, y^. Thus, to maintain these biases at a few centimeters (or 
less than 1 percent) requires that y^ be less than or equal to -50 dB. 

Since -50 dB average RMS sidelobes levels can be achieved with the newer 
pulse compression devices, such as the reflective array compressor (RAC), 
the biases due to range sidelobes should be negligible for the MLE processor. 
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4.0 


split-gate trackers 


The performance of a split-gate tracker is examined in this section 
of the report. The variance of a split-gate tracker was derived previously 
in reference (1). In Section 4.1 the variance of a split-gate tracker is 
rederived to provide a much simpler expression for interpreting tracker 
accuracy. The effect of changing the gate widths and the placement of the 
early gate is examined. A comparison of tracker accuracies for several 
adaptive and nonadaptive split-gate trackers is given for typical SEASAT-A 
parameters . 

Section 4.2 examines the altitude bias caused by range sidelobes in 
the pulse compression system. It is shown that receiver weighting to reduce 
peak sidelobes will not be necessary for SEASAT-A, since peak sidelobes as 
large as -10 dB only produce an altitude bias of about 1 cm. However, the 
far out sidelobe level must be less than -50 dB RMS in order to keep the 
altitude bias less than 6 cm. 

Section 4.3 describes the altitude bias caused by the additive 
transmitter noise burst resulting from gating the TWT on prior to (and keeping 
it on after) transmitting the LFM signal. It is shown that this transmit noise 
burst produces uniform receiver range sidelobes extending over the entire side- 
lobe region. To maintain these sidelobes at -50 dB requires that the ratio 
of total transmit signal energy to total transmit noise energy be at least 
17 dB. 


4.1 Tracker Accuracy 

In this section the performance of a split-gate tracker is derived, 
and the effect of changing the gate widths and the placement of the early gate 
Is examined. To keep the derivations relatively simple, a sampled system will 
be assumed. A split-gate tracker estimates the epoch of a signal by balancing 
the output of an early gate with that of a late gate. Specifically, let: 

t , i = 0, + 1, +2, ... be the times at which the signal is sampled. 
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and 


v(t i ) 

£ 6(t -t)h 
i 


be the square law detected output of 

the receiver at time t., 

1 

be the impulse response of the early 
gate. 


C £ 6(t.-t)g. be the impulse response of the late 

i 1 1 gate 


c 


be a constant which makes the tracker 
unbiased. 


Note that the Kronecker delta functions arise from the assumption of a sampled 
system. 


For any time delay, Tq , between the tracker estimated epoch and the 
true epoch, the error voltage of the split-gate tracker is given by: 

C(T 0 ) = E V(t 1 -T 0 )(h i -Cg i ) (4.1) 

The tracker feedback loops function so that the error voltage is held at 
zero. Thus the tracker estimated epoch is given by the solution of: 


e(T Q ) = 0 (4.2) 

If the loop time constants are long, and the tracking error is 
unbiased, then t q will be nearly zero. Expanding £(Tq) in a Taylor series 
about Tq = 0 and substituting into (4.2) yields: 


T «, -e<o> 

o a&(T) 

dT 


(4.3) 


T=0 


Let E 


^( T ) 

St 


be the mean tracking error slope. Using (4.3) one has: 


T=0 


Var T 


0 


Var 6 (0) 



Y 

E St 
V 

T=0 ) 


(4.4) 
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The tracker will be unbiased if 


E S(0) = 0 


Thus the normalizing constant in equation 4-1 is given by: 


C 


2 V(t.) h. 

£ v(t 4 ) 8l 


where V(t i ) = E V(t i > is 
The mean error 


the mean receiver power output at time t^. 
slope is given by 


(4.5) 


9ES(T) 

St 


h r c Bl > 

T~0 


And the variance of the tracking error signal is given by: 


(4.6) 


Var £(0) = 



V(t.)) (h i -c 


8 i ) 


(4.7) 


where the samples are assumed 
i/j). For Gaussian processes 


to be uncorrelated. (i.e., (V.-V, )-0, 

(exponentially distributed power). 


Var V(t t ) = (e V(t i )) 2 = V 2 ^) (4.8) 

Thus the tracking error variance is given by combining (4.6), (4.7), (4.8) 
and (4.4) to yield: 


Var t 


0 


£v 2 ( ti )(h.-c g.) 2 

V(t 1 )(h.-C g t )) 


(4.9) 


From equation 4.9 it is seen that the tracking error variance depends 
on the shape of the early and late gates and on the characteristics of the mean 
power return, V(t). To simplify the evaluation of equation 4.9, the 
approximation sketched in Figure 4.1 is assumed. 
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Figure 4.1 Sketch of the Average Returned Power Level and the Gate Positioning 



That is , the following simplifications are made 


i) Square - non-overlapping gates are assumed. 

ii) The late gate is positioned on the plateau. 

iii) The early gate is positioned on the leading edge. 

iv) The signal return power is modeled as a linear ramp 
rising to a flat plateau. 

v) Hie receiver noise level is unity. 

If the following definitions are made: 

At = time between independent samples 

a = signal-to-noise ratio 

T = rise time of the leading edge 

T = time duration of the early gate 
E 

T t = time duration of the late gate 

Li 

V = average power return at the center of the early gate 
E 

Np = number of pulses averaged 

0 = = wave height parameter 

h 

CT h = RMS wave height 
O' = .3227 

the mean power return on the ramp may be expressed as 

V(t) = V_ + a(*0t (4.10) 

b 

where t is referenced to the center of the early gate and the constant a has 
been chosen such that the ramp approximates the normal integral in a mean 
square sense. Thus, the mean slope of the ramp 

V(t) = a»3 (4.11) 

and the term in the denominator of equation 4.9 


E VO^Xh^C gi ) = £ V(t t )h. 


(4.12) 
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since V(t^) is zero on the plateau, 
the sum by an integral: 


1 ' Si 


But using equation 4.11 and approximating 



a»pdt 



(4.13) 


Similarily, the numerator of equation 4.9 is given by. 


Z V 2 (t i )(h i -C gj = Z V 2 (t i )h i 2 + Z V 2 (t i )C \ 2 (4.14) 

since the gates are non-overlapping. Then substituting equation 4.10 and 
approximating by an integral yields: 


Ev 2 (t i )h i 2 ~ h 


I 


_E 

2 


(V E +aapt) 4 dt 


JS 

2 


if - 2 2 T 

= At \ E T E + 12 


and since the return is constant on the plateau 


(4-15) 


EV 2 ^)^ 2 .^(l + a) 2 C?I L 


(4.16) 


Finally, the normalizing constant given by 
equation 4.5 


Z V( ti )h. 
ZV(t i )g i 



2 

(V_+aQ'Pt)dt 

£< 

(l+a)T L 


(4.17) 


VT„ 
E E 


(l+a)T 
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Substituting the expressions given by equations 4.13 through 4.17 in equation 
4.9 provides an expression for the tracking error variance on a per pulse basis 


or 


2 

T 

O 



(acrp ) 2 



If the split-gate tracks the kth power point. 


(4.18) 


V-, = 1+ka 
E 


and if N pulses are averaged in the tracking loop 

p 



(4.19) 


(4.20) 


Finally, converting from delay to range, the RMS altitude error (in meters) 



(4.21) 


where At can now be interpreted as a range resolution cell (since the samples 

are assumed to be independent) and T and T the range extent of the early and 

EL 

late gates, respectively. 


An examination of equation 4.21 reveals the following characteristics 
associated with the RMS altitude error of a split-gate tracker. The error 
increases linearly with RMS wave height 0^. It decreases as the number of resolu- 
tion cells in the early and late rate increases and as the tracking point k is 
made smaller. 


Thus, for a non-adaptive split-gate tracker (a tracker with fixed 
early and late gates), the quarter power tracker, k=.25, performs better than 
the half power tracker (k=.5) as shown in Figure 4.2. In this case, the 
maximum width of the early gate is limited by the smallest RMS wave height 
to be encountered. Thus under smooth sea conditions, the early gate can only 
be made one resolution cell wide. Furthermore, with a one resolution cell 
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RMS Estimation Error (cm) 



(meters) 

Figure 4 # 2 Comparison of Tracking Accuracies of Split-Gate Trackers and Adaptive Processors 
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early gate, little is to be gained by marking the late gate much wider 

than 16 resolution cells. That is, increasing the late gate from 16 to 

an infinite number of resolution cells only reduces a by 3%. 

K 

Improved tracking performance can, however, be achieved with 
an adapting split-gate tracker. This is accomplished by estimating the 
RMS wave height and increasing the width of the early gate as the estimated 
wave height increases. The early gate width can be increased in discrete 
steps as suggested by MacArthur (2) or continuously as illustrated in 
Figure 4.2. The performance comparison of the split-gate trackers given 
in Figure 4.2 was based on typical SEASAT-A system parameters. The MLE 
performance Is included in the figure to Indicate the ultimate improvement 
achievable with an optimum tracker. 

4.2 Range Sidelobe Bias 


In a pulse compression system, the major effect of amplitude and 
phase errors is to increase the level of the range sidelobes of the system. 
These sidelobes in general do not significantly affect the altitude tracker 
variance, or the accuracy of estimating wave height. They do, however, 
cause bias errors in the altimeter. In this section, an approximate 
expression is derived for the bias errors due to range sidelobes. The 
expression is then evaluated for typical SEASAT-A parameters. It is shown 
that peak sidelobes as large as -10 dB only produce an altitude bias of 
about 1 cm, thus receiver weighting to reduce the peak sidelobes will not 
be necessary for SEASAT-A. It is also shown, however, that the average far 
out sidelobe level must be less than -50 dB RMS in order to keep the altitude 
bias less than 6.0 cm. 


4-2.1 Theoretical Development 

The mean power versus time output from an altimeter can be written 
as : 


-f 


V (T ) = / CT(t) I P (T - 1 ) I dt 


(4.22) 
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where 


a (t) is the distribution of target cross-section with 

time (i.e., sea surface impulse response) 

p (f) is the normalized radar system point target 

response 

and V(t) is the mean power output. 

The standard split-gate tracker adjusts the epoch so that the early gate 
output exactly balances one-half the late gate output. 

Thus, the tracker finds t such that 
’ o 

V(t o ) = % V(t + A) (4-23) 

where A is the separation between the early and late gate. 

Now, for typical pulse compression systems, the point target response 
consists of a narrow main lobe response plus extended sidelobes as sketched below: 
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(4.24) 


2 

Thus we may decompose the response function |p(t)| as 
2 

Jp(t)| “ + PgL^ 

The mean power output is given by: 

00 

v(r) = J CT (t)p ML (r-t)dt 
“ 00 

+ y^a( t )3 SL ( T_t ) dt (4.25) 

-00 

As long as the wave height is large compared to the system range resolu- 
tion, B mt ( t) acts as essentially a delta function, thus 


V(t) « t c q(t) + V gL (r) 


(4.26) 


where 


X 

c 


v si, (T) 


is the compressed pulse width and 

is the distortion in the response due to 
the presence of the sidelobes. 


That is, 

00 

V S1 <’> ■ / a(t)e sl (T-t)dt 

-CO 


(4.27) 


From (4.26) and (4.23) we have at the estimated epoch: 

% a<V + v sl (t 0 > - % [ t c o(V A) + v SL < V A)] (4.28) 

Now for small bias errors. 


oi t Q ) « cr(0) + d(0)t Q 


(4.29) 
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and 


cr(t 0 + A) ~ ct(A) + a(A) fc 0 


(4.30) 


Where we assume the true epoch is t Q =0. Further, at the true epoch. 


g-(O) — \ a (A) (4.31) 

Applying (4.29), (4.30), and (4.31) to (4.28) yields: 

t 0 |r c (o(°) - 4 <5(4))] - % ^ S1 (t 0 + 4) - < 4 - 32 > 

Now, in general, V CT is a slowly varying function, and we may write 

u 1_J 

approximately : 


V SL (t o> 


V SL< 0) 


(4.33) 


and 

V V 4) ~ Vto) 

So that finally, the bias is given by: 

* V sl ( A> - V (0) 

° T c Ca(0) - ka(L)l 


(4.34) 


(4.35) 


We assume an Impulse response corresponding to a pencil beam antenna as 
given by Barrick (3). That is: 


CT (t) = k $(t/t h ) 



(4.36) 


Where 


$(LCX) 


LCX 



2 

dt 


-00 
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is the normal integral, and 



is the (one sigma) pulse stretching due wave height and 
tj. = .125 (h( 1 +H/a)^ 2 ) 

is the duration of the impulse response for beam width i|j and altitude H, 

B 

where "a" is the earth radius. 

In general, t^ » t^, so that we may write 

. 1 -%(t/<o 2 

‘ /2Wt h e < 4 - 37 > 

And, the bias is given by 



V SL W) - V SL<°>) 


(4.38) 


where we have assumed A s 3t. . 

h 

4.2.2 Evaluation of the Bias 

To evaluate (4.38), we will assume that the sidelobes consist of 
a decaying term plus a uniform level term. The decaying term follows a 
power law, thus : 


where 

and 



T c is the compressed pulse width 

is the uncompressed pulse width 
Y 0 is the (nominal) height of the first sidelobe 


(4.39) 
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is the "far-out" sidelobe level 

n is assumed to be even since it is the decay law 
for the power sidelobes. 


For simplicity, we assume the uncompressed pulse length, » t^» 
so that the error function becomes essentially a step function when compared 
to the sidelobe envelope, thus the impulse response is approximately a pulse 
of duration t^.. 


Substituting (4,36) and (4.39) into (4.27) and making the above 
approximation yields: 


v sl ( t) ■ Y„<i.5) n / (f)" it*/ 1 (t)" at + ^ 

T T C 

c c 




(4.40) 


Where for simplicity we have assumed that > t^ and note that the term 

zero when r = 0 # 

T 


Substituting (4.40) into (4*38) yields: 




Y 0 U.5) 
(n- 1) 


n 




n 





Thus the total error is the sum of the error due to the decaying 
sidelobes, and the error due to the uniform sidelobes. 
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In most cases of interest, ^>50 T c> so that we may neglect terms 

of order T / t T . Secondly, for purposes of specifying a sidelobe budget, we 
c I 

will assume an upper bound roll-off n = 2. Thus, 





t (<- 5 > 


2 




(4.42) 


Table 4.1 gives the altitude bias due to the close-in sidelobes 
versus peak sidelobe for typical SEASAT-A parameters. Table 4.2 gives bias 
due to the far-out sidelobes. 


Table 4.1 


Bias in Centimeters Due 

to Close-In Sidelobes 


t T /x = 1000, 

o 

II 

4^ 


I c 

C 


Peak Sidelobe 

Wave -Height 


<2 Y.) 

°h - - 25 m CT h 

“ 5.0 m 

-10 (dB) 

.07 (cm) 1 

.52 (cm) 

-12 

.04 

.97 

-14 

.03 

.60 

-16 

.02 

.38 

-18 

.012 

.24 

-20 

.008 

.15 


Table 4.2 

Bias in Centimeters Due to Far-Out Sidelobes 


t T /T = 1000, A/t = 44 
I c c 


Average 



Sidelobe Level 

Wave -Height 


°h - - 25 “ 

a, = 5 . 0 m 
h 

-30 ( d B) 

32.0 (cm) 

632 (cm) 

-35 

10.0 

194 

-40 

3.2 

64 

-45 

1.0 

20 

-50 

.3 

6 


1-4-15 




4.3 


Transmitter Noise Bias 


In the recommended SEASAT-A system, the TWT must be turned on 
prior to (and kept on after) transmitting the LFM waveform, so that the 
transient response of the TWT is not superimposed on the desired signal. 
This results in a burst of noise being transmitted before and after the 
signal. This degrades the sidelobe response of the receiver as shown in 
Figure 4.3. To maintain these far out receiver sidelobes at -50 dB (alti- 
tude bias less than 6 cm) requires that the ratio of total transmit signal 
energy to total transmit noise energy be at least 17 dB. 

4.3.1 Analysis 

Let S(t) be the desired transmit signal, and Tlj(t) be a noise 
burst on transmit. In general T1 ^ ( t) will have a time varying power level 
as sketched below. 


vuwm 



Noise 

Burst 
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Thus the transmitted signal may be represented by: 


Z T (t) » S(t) + T)(t) (4.43) 

Let p(t) be the target (distributed) reflection coefficient at range 
cx/2. Then the received signal may be represented by: 


Z R (t) = J* p(t) Z,j,(t-T)dr (4.44) 

T 

If the receiver is not matched, but has receiver characteristic w(u), 
then the (complex) amplitude of the receiver output at time t is given by 

X(t) - f w*(u)Z R (t+u)du (4.45) 

Substituting (4.43) and (4.44) into (4.45) yields. 

X(t) = f p(t) v (t-T)dT + P p(-r) N(t-x)dr (4.46) 

where v (t) is the cross ambiguity function given by (zero doppler) : 

A sw 

X gw (t) - J w*(u)S(t+u)du (4.47) 

and N(t) is the noise process which results from passing T»(t) through the receiver. 
Thus 

N(t) - J* w*(u)Tj(t+u)du (4.48) 
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The mean receiver output is then given by: 


P R (t) - E X*(t)X(t) 

* ]“ g(t) {j Xsw (t-T )| 2 + EjN(t- T ) | 2 | d T 
where a(r) = Ep(T)p*(r) is the mean cross-section of the target with 


If the noise pulse is time-varying white, we may write 


E ('n( t 1 )'fi*(t 2 )) = yt r 


Sin ^(tj^-tj) 

"VW 


where P„(t) is the noise power at time t t 
N 1* 

Thus 


E|N(t)J =; / J w*(u 1 )w(u 2 )ETl(t+u 1 )T|*(t-Hi 2 )du 1 du 2 


U 1 U 2 


= J f P N (t+u 1 ) SinnB N (u r U 2 ) w*(u 1 )w*(u 2 )du 1 


du„ 


U 1 U 2 


"VW 




! 


|„(u)| 2 du 


B 


N 


“ J N^(t+u) fw(u) du 


where is the power spectral density. 


Substituting (4.51) into (4.49) yields 

P R (t) = J* o(t) Us*/ 11 -'!’) | 2 + J* N Q (t-T+u) jw(u) | 2 dul dr 

l n J 


(4.49) 

range. 


(4.50) 


(4.51) 


(4.52) 
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Thus the effect of the noise burst is to additively distort the 
idealized ambiguity function by the term: 

D(t) = J N q (t + u) |w(u)| 2 du (4.53) 

Let P^(t) be the noise power (in the signal band) at time t, 
(relative to the signal). Then N Q (t) = P^(t)/B where B is the signal 
bandwidth. Thus 

D(t) = ^ J P N (t + u) Jw(u) J 2 du (4.54) 

4.3.2 Approximate Evaluation 

To evaluate (4.54) approximately, we assume that the noise burst 
is completely suppressed by the signal, when the signal is present. Further, 
we assume the envelope of the noise burst is square, and we assume that the 
receiver weight is uniform. Thus, the situation is as sketched below: 
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For these assumptions, 


the distortion is given by 


|t| < AT/2 

AT/2 < |t| < T/2 

T/2 s I C 1 < T /2 (4.55) 

1 n 

where AT = T - T 
n 

Referenced to the signal, the distortion term results in the more or less 
uniform sidelobe structure shown below: 


D<t> - 


t Ia 

BT 


AT _N 
2 BT 





Figure 4»3 

Receiver Si delobes Resulting From Transmit Noise Burst 
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Then the maximum sidelobes 


SL 



2E BT 
8 


(4.56) 


where * total noise energy 

E = total signal energy 
8 

As was shown in Section 4.2, far out sidelobes should be kept at less than 
-50 dB to maintain an altitude bias of less than 6 cm at large RMS wave 
heights. Since for the recommended SEASAT-A system 2BTN2000, to achieve a 
-50 dB sidelobe (from equation 4.56) requires that 

E 

^ > 17 dB (4.57) 

e n 

While this requirement should not be difficult to achieve in practice, it 
does point out the fact that care should be taken to minimize the time the 
TWT is gated on prior to (and left on after) transmitting the LFM signal. 
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5.0 


ALTIMETER POWER RETURN 


In this section, the radar equation for a distributed target is 
derived and applied to the particular case of a satellite altimeter. The 
target is modeled as a continuous distribution of independent scattering 
points distributed in the three dimensions (T, v, u) . Here, T is range 
delay, v is the Doppler shift, and u is the sine of the angle to the 
scattering point measured from the antenna boresight. Since the scattering 
points are assumed to be incoherent, the total power is simply the integral 
over delay, Doppler shift, and angle of the differential power received from 
each point. Although the altimeter equation has been derived more rigorously 
elsewhere (e.g., Barrick (3), (4) or Harger (7), this derivation is included 
for completeness, and to define the notation used to obtain an approximate, 
closed-form solution for the shape of the mean power return. 

5 . 1 Theoretical Development 

From the theory of high resolution radar (e.g., Rihaczek (8) 
or Deley (9), the power received from a mls-matched filter can be written 
in terms of the cross-ambiguity function. That is, let the received narrow- 
band signal be represented as: 

S r (t) =/iF , g(t-x) exp[j2lI(f o -v)(t-T)] (5.1) 

where T is the round-trip delay to the scatter, 

v is the Doppler shift produced by the 
scatterer , 

is the carrier frequency, 

g(t) is the complex modulation impressed on 
the carrier. 


and 


is the received signal energy. 



If the receiver is matched to a complex modulation w(t), a Doppler frequency 
0, and a delay, t, then the video power out of the filter at time t is given 
by: 


where 



2 

v(t) = %E r G o B lx gw (t- T > *- v )l 
is the normalized cross-ambiguity function defined by: 


(5.2) 


and 


V T>V) 



g(5) w*(?+T)exp[-j2nvF] d? 



J left )! 2 


dt 


(5-3) 


is the product of receiver gain, times its noise bandwidth. Note that Xg„ 
has been normalized so that 


If the scattering point has differential cross-section, da, and is 
located at angular position u (in sine theta space), then the signal energy 
can be computed from the standard radar equation (Skolnik (10). That is: 


E 


r 


E T GA(u) 

4nR 2 


da 

4I1R 2 


\ GA(u) 
411 


(5.4) 


where 


and 


E t is the transmitted energy, 

G is the peak antenna gain, 

A(u) is the one way loss factor when the target is 
not on boresight, 

R ct/ 2 is the range to the scatterer, 

X is the R.F. wave length 

da is the differential cross-section at the 

scattering point, and is generally a function 
of all three parameters (t, v, u) , 

L is the total system losses (except for 

processing mis-match which is included in 
the cross-ambiguity function). 
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Combining Equations (5.2) and (5.4), gives the differential power 
received from a scattering point located at (T , v , u) : 


dV(t) = 


(tej) g 2 x 2 lg o b 

(4fl) 3 


A 2 (u) 

R 4 


X a „(t-T, *-V)| da (T, V, u) 
gw 


(5.5) 


Integrating this equation over the three variables, and adding the receiver 
noise, yields the total average video power for a distributed target: 


2 

| X (t-T,$-u )| 2 da (t, V, u) 

R gW 

(5.6) 

where k is Bolzman's Constant, 

T is the reference temperature, 
o 

and F„ is the receiver noise figure. 

N 

Equation (5.6) is a fairly general representation of the power received 
from a distributed target. To apply it to the satellite altimeter, one must now 
compute an expression for the differential cross-section da(T, v, u) which is 
appropriate to the altimeter. Hie geometry being considered is shown in 
Figure 5.1. 

Consider a scattering point located at height, h, above the mean sea 
surface. Then, if a°(fl) is the mean cross-section per unit surface area (at 
angle 9 relative to radar) and if P^O 1 / 0 ^ is the probability density function 
of the distribution of scatterers with height, (scaled to RMS wave height o^) 
then the differential cross-section at that point Is: 


V(t) = 


%e t g 2 \ 2 lg o b 

(4II) 3 


/ f f 

U~U T="f°° 


+ %kT F..G B 
z o N o 


da = a°(e) P h (h/a h ) dv (5.7) 

where dV is the differential volume at the point of interest. Now, in the 
spherical coordinate system (R,0,i|f) located at the satellite, the differential 
element of volume is: 
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dV * dR (RdG) (R sin 0 di^ ) 


= sin 0 dR dB di|r (5.8) 

The final step, then is to compute T, v, u and h in terms of R, 

0 and , substitute the results into Equation (5-6) and then substitute 
Equation (5.8) into Equation (5.6) to obtain the final expression. But from 
Figure 5.1, it is seen that 

c 

2V 

v ~ — sin 0 cos i|r 

u = sin 0 (5.9) 

and 

h = ( (H -R cos 6 + a ) 2 + R^sin^ 0 ) ^ - a (5 . 10) 

where V is the satellite velocity, 

H is the satellite altitude, 

o 

and a is the radius of the earth, 

e 

Before making the indicated substitutions, the equations will be 
simplified by making some approximations. That is, it is assumed that the 
beam width of the altimeter is narrow (less than 5°) so that small angle, 
and related consistent approximations can be made. Thus, it is assumed that: 

u ^ sin 0 « 1, 

R = H + AR, 
o 

AR « H 

o 

AR « a (5.11) 

e 

With these approximations, one can write: 
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T ,!H£ + 24B = t + 2A£ 

c c o c 


H 0 H 

h ** -AR + —^ 2 — (l+~) 

e 


( 5 . 12 ) 


da « a°(0) P. (h/a. ) H 2 6 dAR d0 d* 
n n o 


and 


2V0 

v « -j— cos 


Thus, the average video power can be written finally as: 


V(t) 


%E G 2 X 2 I£t°G B 
T O 

(411) 3 H 3 (1+H /a ) 

o v o e -oo 


00 fl/2 
J dAR J 


[H Q (l+H o /a e )e d0] • 


H 0 


V(0)A 2 (0)p, (- & + (1 + H /a )) 

"• \ 0^ 2c ^ o e7 


2n 

/ 


| m 2AR , 2V0 . . 1 2 

1 X gw (t " T o " ~ ’ 0 ‘ T COS ♦>! 


+ %kT F G B (5.13) 

o n o 

where Y(0) is defined by cr°(0) = q°y( 9) and defines the variation of o° (0) 
with angle. 

5.2 Approximate Evaluation 

A closed form solution for the mean power return can now be 
obtained from the triple integral of Equation 5*13. To achieve such a 
solution requires some additional approximations and simplifying assumptions. 
First assume that the cross ambiguity function is essentially constant over 
the doppler variation of the scattering points, then 
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(5.14) 


2tt 

f d tl X gw (t-T o - 0 - COS ♦) |‘ 

•'A 


“H x g „<t-T 0 - 0)| 2 

Further assume that the variation of the cross ambiguity function 
with range delay can be approximated by a Gaussian function, thus 

* “pl-ifr) 1 - /Sr °t T, (ir) <5 ' 15> 

T T 

where T| is the normal density function and <7^ ~ RMS pulse width. 

Now from Barricks (3) model, the wave height density is also 
normal, hence 



‘ 11 % 

(5.16) 

Letting 

H 



H = H (1 + — ) 
0 a e 

(5.17) 

and 

He 2 
x = — 



dx = H0d 9 

(5.18) 


and substituting equation 5.14 through 5.18 into the triple integral of 
equation 5. 13, denoted I^, yields 

«r 2 

Q 

.3/2 


Htt 

8 . 


+CO 


I 3 * (2tt)' 


■. 7 


dxY(x)A (x) / dAR * 


/ 


0 


/ c(t- T o) _ Ar \ 

’ ' K — ■) 


(5.19) 
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where 


C(7 T 

<j R = — - RMS range resolution 


Thus the integral over AR is simply the convolution of two normal densities 
having different means and variances. As may readily be shown, this convolu- 
tion produces a normal density with a mean equal to the sum of the means and a 
variance equal to the sum of the variances. Letting 


C (t-T ) 


y = 


» range measured from the mean sea surface (5.20) 


and performing the indicated convolution of equation 5.19 

2 


HTT 


’ <2 " >3/2 ° r / 


8 


Y(x)A ( 


x)H^ 


dx 


(5.21) 


where 


2 2 2 
a = cr + a, 

e R h 


(5.22) 


Now assuming that the two way antenna pattern function can be 
approximated by a Gaussian function 




(5.23) 


where 


0..,, = the two way 3 dB beamwidth 

DWZ 


And from Barricks (3) model the surface shaping function 

Y(e> r *' ln (\) 


(5.24) 


where 


B g - the 3 dB spread of the slope distribution 
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Then utilizing the change of variable given by equation 5.18, the product of 
the surface shaping function and beam pattern can be written as 

G (x) = Y(x)A 2 (x) = e _8ln 2 HQ 2 (5.25) 

e 

where 

9 = 

e 

is the effective 3 dB beam width. 

The remaining integral of equation 5.21 can now be expressed as 


1 1 
2 + e 2 


-1/Z 


l 9 BW2 


s 


(5.26) 


where 


F(y) * (G * T|) [x] 



(5.27) 


(5.28) 


To evaluate the remaining integral requires completing the square in the 
exponent, let 


%(^) 2 


+ |lX 


% 


2 2 
y -2xy+x 

2 

a 

e 


+ % 



= 1 - , [y 2 -2x(y-a V)+x 2 ] 

2 a 

e 

= [y 2 -(y-a ^J-) 2 + (y-a 2 ji) 2 - 2x(y-a 2 jj,)+x 2 ] 
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1 [y 2 -y 2 +2ycr 4 H 2 ] + [x-(y-a ? Ul 


2.,,2 


2a 


e e 


20 


(yv-%^ 2 ) + % 


x-Cy-C^M')’* 2 


] 


Thus 


F(y) = exp[%(a e M-) -yM- ] 


, l /-» f /*-(7-a,W 

/ exp [ l - ^ - ) 


dx 


Letting 


t = 


(y -a m>)-x 
e 


F (y) = exp[%(a^iJ.) -yM-] 


— 1 - - a p 
a e 
e 


/ 


exp(-%t 2 ) dt 


So that finally 


F (y) = exp[%(a e ti) 2 -yii] - a e ^)J 


(5.29) 


(5.30) 


(5.31) 


where $ is the cumulative normal distribution function. Thus F(y) describes 
the shape of the mean power return as a function of range measured from the 
mean sea surface. 


Substituting equation 5.31 into equation 5.21 and the result into 
equation 5.13 yields the approximate closed form solution for the average 
video power 


V (y) 


%(2Tr) 3/2 E T G 2 \ 2 LG o Ba R 

3 3 Hrt 

(4tt) j H o (1+^) 


e 


• a°(a h )exp[%(a e u.) 



+ %kT F G B 
o n o 


(5.32) 
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where 


10 log ^a° (c^)) = 13,6 - 5 log (5.33) 

is the variation in cross section per unit area of the sea at normal incidence 
as a function of RMS wave height (in meters) as given by Barricks (3) model. 

Figure 5.2 illustrates the shape of the mean power return as given 

by equation 5.32 as a function of RMS wave height for typical SEASAT-A 

parameters. The factor of 4 or 5 variation in power at the so called "plateau” 

of the return can be attributed to the change in cross section as given by 

equation 5.33. If this variation in a°(a n ) is normalized (as with a late gate 

on the plateau) the mean power return of Figure 5.3 is obtained. The droop in 

the plateau region shown in Figure 5.3 is caused by the rather narrow 

effective beam width 0 . A small 8 also causes a distortion in the shape 
e e 

of the leading edge, producing a wave height dependent bias. As shown in Figure 5.4, 
a bias of 19 m occurs at a 5 m RMS wave height due to this effect. A plot of 
distortion and bias as a function of antenna beam width, Figure 5.5, shows that 
an antenna beam width of 5°to 10° would essentially eliminate this problem. 

The bias as a function of antenna beam width for a 5 meter RMS wave height 
is shown in Figure 5.6. While a wider antenna beam width would also reduce the 
effects of antenna pointing errors, the reduction in gain would require a 
corresponding increase in peak transmitter power and/or compression ratio to 
achieve the recommended signal-to-noise ratio. 

The scope of the present program precluded an in depth study of the 
effects of pointing error on the shape of the mean power return (although it 
is presently felt that the effect of pointing error could be included in 
the closed form solution given by equation 5.32). It is recommended that a study 
be made of all the trade offs involved in order to determine an optimum antenna 
beam width . 
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Figure 5.2 Mean Sea Return 
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Figure 5.3 Normalized Mean Sea Return 
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Figure 5.4 
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Figure 5.5 Distortion Due to Antenna Beamwidth in Mean Surface Return 
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Figure 5.6 


Bias Due to Antenna Beamwidth in Mean Surface Return 
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APPENDIX A 
(Part I) 

THEORETICAL DERIVATIONS FOR THE MLE ANALYSIS 


A.O INTRODUCTION 


In Reference (1), a maximum likelihood estimator (MLE) was pro- 
posed which would simultaneously estimate the epoch , wave height, and 
signal-to-noise ratio for a satellite altimeter ♦ In that reference, the 
asymptotic performance, when many pulses are integrated, is analyzed* In 
Reference (2), a technique for implementing the processor using feedback 
loops is proposed, but its performance is not analyzed. In this appendix, 
a somewhat broader class of estimators is considered. That is, those esti- 
mators which are obtained by minimizing a penalty functional which depends 
on both the observed data and on the parameters to be estimated. Maximum 
likelihood estimators (MLE) are special cases of this class of estimators. 

In this appendix, formulas for computing the asymptotic bias and 
variance of the estimators are derived. In addition, an appropriate imple- 
mentation scheme is discussed which is based on feedback loops. Finally, 
formulas for computing the mean response, and the variance of the loop 
estimators are derived. The general formulas are then evaluated for the 
satellite altimeter case. The resulting formulas provide the basis for 
the results discussed in Section 3.0, 

A. 1 Problem Definition 

It is assumed that for each pulse (or data interval), a vector 
of K data values are measured. Thus, let: 

be a K dimensional vector of measurements 
taken on the p U i se4 

For the altimeter problem, the vector represents the video outputs from 
K range cells. For a given set of parameters, represented by the M dimen- 
sional vector 0, a penalty function is defined, call it: 

q ( v 1 /e> 


A-l 



Then, for an entire sequence of N pulses, the total penalty is 


N 

£ Q(v./0) 

i - 1 

A 

The class of estimators to be considered are those estimates 9 which mini- 
mize the total penalty, i.e., 

N A N 

Z Q(V,/e) - min Z Q(V./0) (A-l) 

i = 1 0 i - 1 

If one lets y T = (S/d0^, , 3/d0^) be the gradient 

with respect to 0, then (assuming existence of the derivatives) a necessary 
condition on the estimator is: 

N 

I CvQ(V./B)] = 0 (A-2) 

1-1 

Before concluding this section, it should be noted that to derive 
a MLE estimate from this formulation, one lets the penalty function be the 
negative logarithm of the likelihood function. That is: 

QO^/B) = -ln(p(V./Q)) 

where P(V./8) is the joint probability density function 

(likelihood function) of the observed values 
given the parameter vector 9. 

Note that minimizing Q is equivalent to maximizing P. Note also, to derive 
a least mean square estimate, one lets 

Q(v./e) = [v. - v(e)] T w[v i - v(e)] 

where V(Q) is a model for the mean value which depends on the 

parameters to be estimated, 

and W is a normalizing matrix. (Note for W = I, one gets 

the usual unweighted least square formulation.) 

A. 2 Asymptotic Bias and Variance 

The asymptotic bias and variance can be estimated by expanding 
(A-2) in a Taylor Series about the true parameter values. Thus, Equation 
(A-2) becomes: 


C. 






A-2 



0 


N N T 

= 2 V Q (V / 9) + E [ V V Q<V./0)3 (9 - 9) 

1=1 1 i= 1 


+ higher order terms. 


Asymptotically, two assumptions are made: 

A 

a) 0 is close enough to 9 so that only the first term 
in the Taylor Series expansion needs to be kept, 

b) Enough independent pulses are included so that the 
second summation may be replaced by its mean value. 
Thus: 



If one lets: 


E [ 7 V T Q(V. /9) 3 » N E V v T Q(V./0) 
i = 1 1 1 


R( 0) = E V / QO^/0) 


and 

e.(9) - v QCV./9) 

Then Equation (A-3) becomes: 

N 

N R(0 ) (6 - 0) = - E e. (9) 
i* 1 1 


(A-3) 


(A- 4) 


(A-5) 


or 



(A-6) 


From (A-6) the asymptotic bias and cross correlation matrix of 0 
is given by: 


E(0 - 6 ) = -R - 1 ( 0 ) Ee(0) 


and 


E (6 - 0) (9 - 0) T = N _1 R _ 1 (0)[Ee(0)e T (9)]R' 1 (0) 


(A- 7) 


(A- 8 ) 
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A. 3 Special Case: MLE 

For the MLE, Q has the special form: 

Q (V / 0) = -In P(V/0) (A-9) 

And thus: 

e(6) - V Q(V/e>- - V (p(V/0)) (A- 10) 

The bias is computed by: 

E e(0) - [yP(V/6>] P(V/6) dV 

= -J v P(V/9) dV 
= -V J* P (V/0) dV 

--71 = 0 (A- 11) 

Thus, by (A-7) the MLE is asymptotically unbiased. It should be 
noted that the above manipulations are formal, and can be justified only 
for regular distributions. 

The information matrix, R(0), is computed by: 

R(0) - E V V T Q(V/9) 

= -E v [P(V/0) T 1 V T P(V/0) 

= E (p(V/ 0))' 2 (v P(V/0)) (v T P(V/0)) 

- E P(V/G) -1 V V T P(V/0) (A- 12) 

Now recall, (Equation (A- 10)): 

€(0) = (pCV/0))' 1 V P(V/ 0) (A- 10) 

Thus: 

R(0) = E e(0) e T (0) - E P<v/0) _1 V V T P(V/0) (A- 13) 

But 

E P(V/0) _1 V V T P(V/0) = J PCV/6)" 1 (v V T P(V/0)) P(V/0) dV 

- 7 V T J P(V/0) dV = V V T 1 = 0 (A- 14) 
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Again, the above manipulations are only justified for regular 
distributions. From (A-13) and (A-14) one has finally: 

R(0) = E e(0) e T (0) (A-15) 

In summary, the asymptotic bias and covariance of the MLE estimates 

are: 

E (0 - 0) = 0 (A- 16) 

and 

E(0- 0) (0 - 0) T = N" 1 [E e(e)e T (0)]" 1 (A- 17) 

where 

e(0) = -V In P(V/0) (A- 18) 

Equation (A- 17) is equal to the Cramer-Rao lower bound, thus, as 
is well known, the MLE is asymptotically efficient, 

A. 4 A Feedback Loop Implementation 

Equation (A- 2) stated that a necessary condition to obtain the 
desired estimator is to find 0 so that the average error signal is zero. 

That is, 0 must satisfy: 

N 

E €. (0) = 0 (A- 2) 

i = 1 X 

One technique for solving this equation is to let the error sig- 
nal be an input to an integrator 

eft = -k E e.(^(t>) (A- 19) 

dt i = 1 1 

If the loops are stable, then at steady state, d'B/dt = 0, and 
the steady state value of 0 must satisfy Equation (A- 2). Hence, at steady 
state *0 converges to the desired estimate 0 which minimizes the total 
penalty. 

The loop equation, (A- 19), can be shown to always converge by 
the following argument. By differentiating the penalty function, one gets: 

d ^ ~ ^ t ® t ~ rfa 

JZ E Q(V./0) = E 9 Q(V,/0)~= E C, (0)~ (A-20) 

i = 1 i = 1 i = 1 
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But, by Equation (A-19) this results in: 


. N N 

E Q(V /e) = -k [ r 
dt i = 1 1 i - 1 



N 


('e)][ e 

i - i 



* o 


(A-21) 


Thus, the total penalty is decreasing with time and must, there 
fore, converge to a local minimum. 


The loop convergence time may be computed approximately by using 
the asymptotic linearization: 


N N 

z 6,(0) » E e. (9) + N R(0) (6 - 0) (A- 22) 

i-1 i = 1 

A 

But, by using Equation (A-6), the first term may be written in terms of 0, 
thus: 


And one has: 


N 

E 

i - 1 


e i ('0) = N R(6) [- (0 - 0) + 'S'- 6] 

= N R(0)(0- 0) 


a 

dt 


-kN R(0)(fS- 0) 


(A- 23) 


(A- 24) 


The solution to Equation (A-24) is: 

t 

0(t) = e" kN R<:9)t: ? o + e -kN R(0)t y* e kN R(0)T (kN R(0)e(T)) dT (A-25) 

o 

where is the initial value for the estimate 0. Note in the above expres- 
sion, 0 may vary during the loop settling time, but it is assumed that 0 is 
not varying* By assumption 0(t) represents the minimum penalty estimate for 
N puises centered at time t* For simplicity, we assume that the estimates 
are only computed at discrete times. Thus: 

0(t) = pff E 6 (t - jN/PRF) (A-26) 
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where 

and 


PRF is the pulse repetition frequency 
6 is the Dirac delta function 

A 

0. is the minimum penalty estimate computed for the pulses 

J occurring from time (j - %)N/PRF to time (j + %)N/PRF. 

From (A-25) and (A-26) one has that the mean of *9(t) is given by: 


E *(t> - e-“ R <® C e o + jfj E R(e >< t ’ J N/PRr > kN R(6) E(5) (A-27) 

J<t2 r 

Assuming that the loop time constant is large compared to N/PRF, then 

E'S(t) „ e ’ kN R(e)t V + e- kH *<»' / s kH R(0 > T kN E( 6 ) dT(Ee) 

o •/ 


= E § + e" kN R(e)t (^ - E 0) 

Thus, the behavior of the mean of 'ft is determined by the eigen- 
values of the matrix: 


(A- 28) 


R(0) = E \7V Q(V/0) 

In particular, if all the eigenvalues of R( 0) are positive, then: 


(A-29) 


lim -kN R( 0) t _ 
t -*» 6 “ U 


and 


E 6 < t ) n rr E e 

t -4 oo 

Further, the convergence time of the loops is determined by the minimum 
eigenvalue of R(8), thus the maximum loop time constant is given by: 

T c = (kN Xmin) (A-30) 

where Xmin is the minimum eigenvalue of R(0). 

The variance of the estimate 'e is computed as follows. From 
(A-25) and (A-28), one has: 

U<t) - E Set) - e-“ R <»‘ / e kN R(9)T kN R(8) (e(T) - E 5 ) dr (A-31) 


A- 7 



Thus: 


Cov ('lift)) - E ( ?<t) - E 'ij(t)) (?(t) - E ?(t)) 1 

.// R(6) (t - Tj) ( kNR(e) ) c<)V (- e(Tl) , e(T 2 )) 

o o 

• (kN R(0)) e" kN R(0)(t- T 2> dT 1 dT 2 (A-32) 
Now, consider Cov ^0(T^), 6 (t 2 ))> from (A-26) one has: 

Cov (e<T.), 0(T )) = fefj) 2 £ Cov [0 , 9 m ] 6 (Tj-tN/PRF) 6 (Tj - mH/PRF) 

't-.m (A-33) 


But, assuming independent pulses, one has: 

jv G 

0 j ± l 


Cov 0. 0. - Cov 0 j = -t 

J ■v 


(A- 34) 


Thus 


Cov 9(T^) 0(T2) = £ Cov(0j 0^) 6 (t^ - jN/PRF) 6 (t 2 “ jN/PRF) 

j ,-L 

= Cov 0 S 6 (T 1 - jN/PRF) 6 (t 2 - jN/PRF) 

Substituting the above expression into Equation (A-32) gives: 


(A-35) 


Cov^G (t) = C^) 2 E e 

J N 


■M R<0) (t - jS/PRT) ( k „ R(0) ) Cov §( mr(6 )) 


. e -kN R(0)(t- jN/PRF) (A _ 36) 


Now if it is assumed that the loop time constant is large compared to N/PRF, 
then 

Cov (So:)) e-“ R(9)(t - T> (kH w cov "0 (kN R) e' RN R ‘®<t-T) dT (A . 37) 

o 

The above equation can be evaluated by expanding exp [-kN R(0)(t-T)] 
in terms of the eigenvalues and eigenvectors of R. Thus: 


e -kR R < 9 > < C - T > M R(0) = 1 e-“ V' - x> T 


1=1 


h *1 *1 


(A-38) 
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where 

and 




is the eigenvalue of R(0) 

is the corresponding eigenvector. 


Thus: 


c s/ 

£,m o 


*V dT \ % tcov 0] cp m cp m T 


kN“ i _ e “kN(>^ + \) t j * T 

E C“ IT" (\V ]< et^ [Cov 6:1 fm 


£ 

PRF 


f,,m % 


W 


(A-39) 


For the special case of the MLE, 

Cov 0 = R' 1 (9) = N _1 2 X.” 1 cp cp T 

j J J J 


Thus (A-39) becomes: 


Cov fl)(t) ] = 2(1- e _2kN C ) cp^ cp/ 


l 


kN /T -2kN R(0)t. 
2PRT (I ‘ 6 > 


(A-40) 


Further, as t -» o° 

Cov I (A-41) 

Equation (A-41) shows that the steady- state variance of ^ can be 
made small by reducing the feed-back loop gain. However, Equation (A- 30) 
shows that this increases the settling time of the loops. Thus the reduced 
variance is achieved by increasing the loop averaging time. 

A. 5 Evaluation of the Satellite Altimetry Case 

In the previous sections, the asymptotic variance of a MLE was 
shown to depend only on the information matrix: 

R(6) = E e(0) e T ( 0) (A- 15) 

where 


e( 9 ) = v Q(v/ 9 ) 


(A- 18) 



and 


Q(V/0) = -In P(V/0) 

where P(V/0) is the joint likelihood function of the vector of observed 
values, V, given the parameter values 0. In this section, the information 
matrix will be evaluated for the satellite altimetry case. 

For this section, it is assumed that the vector 


V - {Vj}j =1 (A-42) 

represents the sampled video output of the altimeter. It is assumed that 
the sampled values occur at least one range resolution cell apart, thus 
the sampled values are statistically independent. Further, square law 
detection is assumed so that the video samples follow an exponential dis- 
tribution. With these assumptions, the likelihood function can be written: 


P(V/0) = 




(A-43) 


where 


V (0> = E (V /0) (A- 44) 

Thus, the likelihood function depends only on the mean power return as a 
function of range. 


Applying (A-43) to (A-18), one obtains the penalty function 

K 

Q(V/0) » E In V. + V./V (A-45) 

j = 1 J 2 

Thus: 

e(0) - 7 Q(V/0) 


= E (V - V.) (V.)” 2 V V. (A-46) 

j = 1 J J J J 


And the information matrix is given by: 
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R(0) = E e(0) € T (e) 


K K 
= E E E 
J-U-l 


(V - V.)(V 4 )" 2 VV, V T V 
J J J J ^ 



- V 


- (Vj -1 V VjHTj' 1 V Vj) 1 


K 

= s 

J-l 


B. B. 
J J 


T 


Where the vector B . 

J 


is given by: 


B. = V. -1 V V. 
1 J j 


(A-47) 


(A-48) 


Note that the summation in Equation (A-47) is over range. In 
general, the range resolution cell is small, therefore the summation can 
be approximated by an integral so that: 


-1 / rp 

R(6) - J B(t) B L (t) dt (A-49) 

T 

min 

where the limits of integration in (A-49) are from the minimum to the max- 
imum range sampled, and is treated as a function of time delay, B(t). 

In (A-49) At is the time in seconds between adjacent range samples. Note 
that the dependence on the parameters is contained entirely in V. 

To evaluate (A-49) further, an expression for V(t/0) is required 
(i.e., the mean power return versus time given the parameter values). In 
Section 5.0 a fairly detailed derivation of the mean power returned from a 
distributed target is given, and a reasonably good approximation for the 
altimeter case is derived. The resulting expression depends on compressed 
pulse length, RMS wave height, antenna beamwidth, and the RMS surface slopes. 
However, from the results in that section, it can be shown that when the 
compressed pulse length is small compared to the wave height, and when the 
antenna beamwidth is broad (>5°), the mean power return can be written 
approximately as: 

V(t/0) = a $ (p(t-T o )) + 1 (A-50) 


A-ll 



where the significant parameters to be estimated are: 


with 


and 


0 = (a, T q , P) 

a = signaL-to-noise ratio 

T = epoch 

o 

B = -r^— is a measure of the RMS wave height. 


In Equation (A-50) $(*) is the unit normal cumulative density 
function. Further, for convenience, the known receiver noise power is 
assumed to be unity. 


with 


With these assumptions, 
B(t) 


one can write: 



B a (t) = a' 1 


f (P(t - T q )) 

$ (P(t - t q )) + 1/a 


B (t) = ft -tO(t-T 0 )) \ 

o \$(P(t - T Q ) + 1/a/ 

/ p(t- T ) $ (P(t- T )\ 

V c) - p \ «(P(t-T 0 )) — TT-) 


(A-51) 


(A- 52) 


If the above expressions are substituted directly into Equation 
(A-50), the resulting integrals are difficult to evaluate* Therefore, to 
simplify the computations, the cumulative normal is approximated by three 
straight line segments. That is: 

$ (u) ^ 0 u < -^a 

au + % | u | <.k<x - 1 (A-53) 

« l u > %a -1 
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where the constant, a = .3227, is chosen to make the straight line segments 
fit the cumulative normal in a least squares sense. 

Using the above approximation, and making the change of variable: 

Z = a p(t-T o ) + % + 1/a (A- 54) 

one finds that the components of the vector B(t) can be written in terms of 
Z as follows: 

Since $ (p(t - t q )) = a p (t - t q ) + % , j 0(t - t q ) | £ %a -1 

= Z - 1/a , 1/a £ Z <; 1 + 1/a 

=1 , Z > 1 + 1/a (A-55) 

and $ |p(t - rjj = a , |P(t-T Q )| <: %a _1 

1/a s Z fil+ 1/a 

= 0 , otherwise (A-56) 

Then B (t) = a~^ [— - , 1/a £ Z £ 1 + 1/a 

a tL 

= a" 1 (1 - 1/a Z) , 1/a £ Z <: 1 + 1/a 

= a -1 (a/a+1) , Z > 1 + 1/a (A-57) 

b t (z) - 0 (- 0 /Z) , l/a s z s 1 + l/a 

O 

B p (Z) = P~ 1 ( " - ^z 1/a ) = (1 - . l/a « Z S 1 + l/« 

Substituting the above equations into (A- 50) and integrating over 
Z gives the results: 

R(e) = SJ& D c D (A - 53) 

where D is the diagonal matrix. 

D = / a -1 0 0 \ 

I 0 P 0 J 

' o 0 f 1 / (A- 59) 
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and C is given by: 


where 


. A 


\ 


aa aT 


aT 


v 


T T 
o o 


V 


J ap 


V 

W. 


Caa 


-/ 


1 + 1/3 /* apT max +% ' 1/3 

dZ+ J 


(1 


1/a 


1 + 1/a 


<TTI > 2 « 


- 1 - | log (l + a) +T 1- + (cxPT^-i)^) 2 


1+1/a 
ar 

o 1/a 




= -a [log (l + a) - 


l + a 


1+1/a 
aP 1/a 


4 


2aZ 


i a + 4 a /i , \ , a + 2 

1 " 2a log (1 + a > + 2(a + 1) 


.1 + 1/a 


C T T = a - 
o o 1/a 


/ i - 


2 2 
q a 

(l + a) 


V 


-4 


1 + 1/a 


a n a + 2, 

Z (1 - ~2aZ ) dZ 


-a [log (X + .) 


(A-60) 


(A-61) 
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.1+ 1/a 


'PP 




1/a 


a+ 2 /, , v , (a + 2) 

= 1 ’ — l0g (1 + a) + 4(1 + a) 

Since the variance of the MLE is given by N ^R ^(0), it would be 
useful to have a general formula for C ^ (from which R ^(0) could be com- 
puted) . Although it is feasible to write down the inverse of a 3 by 3 
matrix, the resulting formulas are unmanageably cumbersome. Therefore, we 
proceed differently. 

Let 


aa 


- 1 log (! + «>+ T^+%^)- 


(A-62) 


Then if 


one has 


C' = 

/ c ’ 

/ aa 

c 

aT o 

C ap \ 



1 CaT 
i ° 

C T T 
0 0 

V j 



\ c *p 

V 

c pp/ 

(A-63) 

c = 

C* + d y 

/1\ (1 0 

0) 

(A-64) 


0 

, 0 j 


where 


d ■ • «<T fi > 2 


(A-65) 


The advantage of writing C in this form is that C' depends only on signal - 
to-noise ratio; therefore, (C') -1 can be easily tabulated. Further, using 
the Woodbury identity, C - ^ can be computed from (C') \ That is, let 

F = (C')" 1 


Then: 




= 


1 + dF 


11 


F F 

11 lj 


(A-66) 


where C 1 ^ is the (i,j) element of C ", Table A- 1 tabulates versus 
signal-to-noise ratio from -10 dB to 30 dB in 5 dB steps. 
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TABLE A- I 


F Matrix vs S/N 


S/N (dB) 

F n 

F 12 

F 13 

F 22 

F 23 

F 

33 

-10 

242.0 

374.96 

-242.0 

1640.5 

-440.0 

1563.6 

- 5 

34.649 

53.686 

- 34.649 

212.78 

- 76.271 

194.19 

0 

o 

• 

00 

12.395 

- 8.0 

41.587 

- 21.408 

33.579 

5 : 

3.4649 

5.3686 

- 3.4649 

15.416 

- 9.8626 

9.978 

10 

2.42 

3.7496 

- 2.42 

9.8201 

- 6.5847 

5.141 

15 

2.1285 

3.2979 

- 2.1285 

8.2012 

- 5.4429 

3.7847 

20 

2.0402 

3.1611 

- 2.0402 

7.6221 

- 4.9907 

3.3143 

25 

2.0127 

3.1185 

- 2.0127 

7.3857 

- 4.7978 

3.1303 
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Example 


Suppose it is required to compute the variance in the 
mate when N = 1500, » .5m and — ~ = 23 m , O n = 5a and , 

Now from Equation (A-58) , 

R -1 (0) = apAt d" 1 c' 1 d _1 

where D ^ = diag (a, p \ P) is a diagonal matrix. 


Now, 




■( 


.3227 cT 


max 


2 cl 


n 



= .4003 


From Table A- I, 

F u = 2.42, F 12 = 3.75, F 22 = 9.82 


Thus 


T T 
, o o 


6.961 


and 


0% — - [N 

Vo 


1 R~ 1 ( 9) 3 t t - (apAt) n 

o o 


■1 rr2 


T T 

. o o 


or 


£ 

2 



«<■¥> ^ 


N 


-1 


T T 
, o o 


- 6 . 1 cm 

A. 6 Bias in the MLE Due to Range Sidelobes 

In the satellite altimeter case, the MLE minimizes the 
functional: 


K 

Q(V/6) = E^ [Vj/VjO)] + In V^e) 


epoch esti- 
= 10 dB. 

(A-67) 

(A-68) 


(A-69) 

penalty 

(A-45) 
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where the summation is over range. In (A-45) , V^(6) is a model of the mean 
power returned versus range. The model, of course, depends on 0 the vector 
of parameters to be estimated. The resulting estimate of 0 is a true MLE 
and will be unbiased only when the mean power model is complete and includes 
all significant parameters. In particular, the mean return model derived in 
Section 5 .0 does not include the effect of range sidelobes. Therefore, when 
range sidelobes are significant, and a model such as proposed in Section 5-0 
is used, the resulting estimate will not be a MLE, and will, in general, have 
biases. However, the estimate is still minimizing the functional (A-45), 
therefore the previously developed theory may be applied to compute the bias. 

Recall from (A-7), the bias is given by: 

E(0-0) = -R _1 (e) E e(6) (A-7) 

where 

R(6) = E W T Q(V/0) (A-4) 

and 

E 6(0) = E V Q(V/0) (A- 5) 


To apply (A-7), (A-4), and (A-5) to (A-45) let 

rv; 

E V. = V. 

3 3 

be the true returned power including the effects of range sidelobes. 
Then, 

K _ _ _ 2 - 

V Q(V/0) = E (V, -V.) V. V V, 

j = ^ J J J J 


(A-70) 


and 


K 


E e(0) = E V Q(V/9) = E (V. -V ) v." 2 V V 

j = 1 J 3 3 3 


(A- 71) 

Differentiating (A-71) again, and taking expected values, one 


obtains: 
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R(e> = e v7 Q(v/e) 


= - e (v " 2 )(v v )(V v ) T 
j - 1 3 J J 

+ E (V. -v.){v .' 2 W T V. - 2V " 3 VV y T V } (A-72) 

j _ ^ J J J J J J J 

If the average side lobe level is small, then one may assume: 

r*t — 

V . - V 

■ i ■ -i « 1 (A- 73) 

V. 

J 

Thus the second summation in Equation (A-72) is small compared to the first 
summation, and one has approximately: 

R(0) w - E B. B T (A- 74) 

j = l J J 

where 

B. = V.' 1 V V. (A- 75) 

J J J 

Note that this is exactly the negative of the matrix which was 
evaluated in the last section* Thus, to evaluate the bias due to range 
sidelobes, one only needs evaluate E e(0) as given by Equation (A-71), 
and then evaluate Equation (A- 7) using the matrix derived in the last 
section* 

To evaluate the mean error signal as given by Equation (A-7), an 
expression for V-V is required. This can be derived by including a uni- 
form sidelobe level in the ambiguity function used in Section 5.0. That is, 
in Section 5*0 the ambiguity function was approximated by: 

!x 2 Ct x > | ~ e‘ %(t l /a r> (5.15) 
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N> | O 


If a uniform sidelobe level is added to the above expression, 
then one has: 





(A- 76) 


where T is the uncompressed pulse length, and y. is the average sidelobe 
u 1 

level. In the approximation (A-76), only the far-out, uniform, range side- 

lobes are considered. This is felt to be a reasonable approach since the 

results in Section 4.2, on the split-gate tracker bias, indicate that the 

far-out sidelobes are the most important ones to consider. 

If the approximation (A-76) is substituted into (5.13) > and all 
the various variable changes required to get to Equation (5.19) are made, 
and if the main lobe contribution is subtracted, one obtains: 

_ %E G 2 X 2 L o° G o B 

V(t) - V (t) = — ^ • I (A- 77) 

(4tt) J H J (l + H /a ) J 


where is the integral: 


Si! 

f 8 

1 3 = 2 TTYi J dx Y (x 
J o 


) A 2 (x) / 


c/2(t - T + T ) 
r O u 


c/2(t - T - T ) 
o u 


\ < a 

V h u h 


where all the various symbols are defined in Section 5.0. Assuming that 
» 0 then the second integral above can be evaluated approximately 

as: 

~(t - t +T ) 

/ 2 o u 

dAR \ (- 2 “— ) « 1 ; | (t - T - T ) s x S | (t - T + T ) (A- 79 ) 

p. a. 2 o u 2 o u 

|(t-T -T ) h h 

Z O U 

Thus the integral 1 ^ becomes: 

x ^ 
max 

i 3 - ( 2 tt) y 1 1 dx y(x) A 2 (x) (A-80) 

x . 
mxn 
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where: 


X = Max [0, c/2(t-T -T )] 
min o u 

X = Max [0, c/2(t-T+T ) ] 
max o u 


By Equation (5.25) y(x) A (x) can be written: 

Y(x) A 2 ( x ) = exp [-(8 In 2/H0 g 2 )x 3 


Thus becomes : 


2tt Y, H0 ‘ 
1 e 

8 In 2 


I 1 ' exp |-TTT (! <t - T o + V)H <A - 81) 


where it has been assumed that only those times for which x . =0 and 

mm 

x ^0 are of interest (i.e,, |t-T I <T ). 
max v * 1 o l u 

Substituting (A-81) into (A-77) yields: 


G 2 X 2 Lo° G B H0 2 

v(t) - v(t) = = — = 2 — (2tt ) Yl -£T-t 

(4n) 3 H Q 3 d + H o /a e ) 1 81n2 


(l - exp | 


8 In 2 
„„ 2 


Now, from Equation (5.32), the signal- to-noi se ratio is given by: 


% E t G 2 x 2 L CT° G o B (2tt) 3/2 G 


(4rr) 3 H 3 (1 + H /a )(%kT F G B) 
^ o o e o n o ' 


(A-83) 


Substituting into (A-82) yields: 


V(t) - V(t) = 


a H0 

e 

/2 Fo Y 1 8 In 2 
R 


h l 1 '^ wll <A - 84) 

rib 

e 


where (as before) it is assumed that the noise power has been normalized to 

unity (i*e. * %kT F G B = 1)* In most cases of interest, it may be 
o n o 

assumed that: 


c<t~ T o) 


e df -1 
8 In 2 * * * 


(A- 85) 
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That is, the range over which the return is sampled is much smaller than 
the effective depth of the antenna beam. Here, the depth of the beam is 
defined as the range change on the earth's surface between the center and 
the edge of the hntenna beam. In the above expression, pT* is the depth 
of the beam divided by In 2. 


With this assumption, the difference in the return due to side- 
lobes can be written: 


V (t) - V(t) = 




-( 

Vr ' 


1 - e 


c(t - T + T ) 
o u \ 

-1 ) 


2H 


[l - exp (- cT u /2 p J ' 1 )+exp(-cT u /2^“ 1 )[c(t- t o )/2 (J l' 1 ] (A-86) 

9 m-r /-t ^ 


ay, M- 

\/2tt a 


Making the change of variable, Z = ap(t - t ) + %+ 1/a, and doing some 
manipulation yields finally: 


V - V m a g x [g 2 + z] 


(A-87) 


where 


\ °h 

2tt a a 

R 



T ) 

u 


(A-88) 


g 


2 



(A-89) 


Substituting this expression into Equation (A-71) gives the following 
expression for the mean error signal: 


K 

e e(0) = s J _ 2 
J-l 




(A-90) 




S 1 

apA t 


/ +1/a g 

J ( 1+ T )B(z)dz+ 7TT 

1/ a 



apT 


max 


+ k+ 1/a 

(g 2 + Z) B (Z) dZ 


1+1/a 
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where B(Z) was defined by Equation (A- 75). The above integrals can be 
evaluated in closed form; the manipulations involved are fairly tedious, 
but straightforward, and will not be repeated here. The results are: 


e a = 


Si 


ap A ta 


If 


l + <8 2 --) 


log(l + a) - 


ial 

1 + a 


k l + a ; 


s 2 (apT max- %) 


+ %( 


a ) 2 (aPT +h+-) 2 
K max a / 


1 + a 


- % 




(1 + a) + 


S 2 3 
1 + a 


e 


P 


ap 2 At 



^ir) log(l + a) - 


g 2 (a+ 2)a 
2(a+l) 


(A-91) 


These equations appear to involve too many parameters to easily 
tabulate or graph them. Therefore, they were programmed, and the bias 
errors were computed for a typical case of interest* For this table, the 
satellite height was 725 km, the (one-way) beamwidth was 5°, and the side- 
lobe level was -50 dB* 

The 5° beamwidth is larger than what is proposed for SEASAT; 
however, the broader beam represents a "worst 11 case in that biases due to 
range sidelobes increase as the antenna beamwidth is increased* Further, 
the results in Section 5 . 0 indicate that the approximation (A-50) is not 
valid since it does not include "beamwidth roll-off". It should be noted 
that in a three loop tracker, which does not compensate for beamwidth 
roll-off, bias errors will develop* These biases are likely to be more 
significant than those due to range sidelobes* Unfortunately, funds on 
this contract ran out before this effect could be evaluated. 

As stated earlier. Table A- 2 shows the effect of range sidelobe 
errors. A surprising feature of the table (in view of the complexity of 
the equations) is that the results are independent of signal-to-noise 
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ratio, and that the epoch and wave height errors scale in proportion to 
RMS wave height. These results have not been verified analytically, but 
if they do hold generally, the biases have the simple form*: 


2 a h 

AT ^ 880 ( -) y- 

O c L 

A a h « 600 a h y x 
A a sa 600 a y^ 


* It can be shown analytically that the biases are proportional to 
average sidelobe level, y . 
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Table A-2 


MLE Altimeter Bias Due to Range Sidelobes 

Average Sidelobe Level - -50 dB 
Antenna Beamwidth = 5° (one way) 

Satellite Altitude = 725 km 


Epoch Error (cm) 


IdZB&nilH 

iwmmsm 

10 

20 

Sig. Waveheight (m) 

0 

-i.i 

-2.2 

-4.5 


5 

-i.i 

-2.2 

-4.5 


10 

-i.i 

-2.2 

-4.5 


20 

-i.i 

-2.2 

-4.5 




Waveheight Error 

( cn 0 


l!nwl351H 

5 

10 

20 


0 

.75 

1.5 

3.0 


i 5 

.75 

1.5 

3.0 


10 

.75 

1.5 

3.0 


20 

.75 

1.5 

3.0 




S/N Error (dB) 

* 


EHn&Knm 

5 

10 

20 

Sig. Waveheight (m) 

0 

.03 

.03 

.03 


5 

.03 

.03 

.03 


10 

.03 

.03 

.03 


20 

.03 

.03 

.03 



* The quantity in the above table is 10 log (1 + ASN/SN) 
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Part II 

OCEAN WAVE SPECTRA 

1.0 INTRODUCTION AND SUMMARY 

In Part I, the design of a high resolution radar for altimetry and 
wave height estimation was discussed. The feasibility of utilizing this same 
radar for the measurement of ocean wave spectra as well is examined in this 
part of the report. A preliminary design analysis shows the resulting system 
parameters for the ocean wave spectrometer to correspond very closely to those 
of the radar altimeter design. 

Section 2.0 examines the implementation and design of a high resolu- 
tion radar for the measurement of directional ocean wave spectra. Commonality 
of the implementation with that of the radar altimeter leads to the recommenda- 
tion of a short pulse wide-band radar technique instead of the mathematically 
equivalent two frequency correlation technique. Both techniques perform a 
spectral analysis of the amplitude modulation (caused by ocean waves) on the 
return signal at a given analysis direction. A preliminary analysis and design 
shows that the required modifications to the radar altimeter design which will 
allow the system to measure wave spectra are that the transmitter output is 
switched to a 2 meter antenna which is steered 20° away from nadir and scanned 
conically with a period of 5 sec. On receive the full deramp processor is replaced 
by a surface wave pulse compressor and detector followed by a band of 11 one-third 
(1/3) octave filters to perform the spectral analysis. The filter outputs are 
then detected and integrated for 26 p>sec per pulse. While this preliminary 
design meets all the specified performance requirements, it is based primarily 
on compatibility with the radar altimeter and probably does not represent an 
optimum configuration. A more detailed study of the trade-offs in achieving 
an optimum configuration is recommended. 

Section 3.0 provides the preliminary analysis of the directional wave 
spectrometer upon which much of the system design is based. The analysis should 
be considered a quick "first cut 11 at the problem because of the many simplifying 
assumptions required in a study of this size. Nevertheless, the overall 
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approach, considerations, and conclusions are believed to be correct. It is 
shown that the input to the spectrum analyzer consists of two noise processes, 
one being the desired modulation of a ° due to the gravity wave spectrum of the 
ocean surface and the other due to the fact that the modulation is riding on a 
noise process whose mean is ct°. Standard statistical theory for the detection 
of a Guassian signal in Gaussian noise is used to develop approximate expressions 
for the signal-to-noise ratio necessary for reliable detection of the modulation 
and resulting design constraint. It is shown that the problem of measuring ocean 
wave length is analogous to the beam limited satellite altimeter problem Results 
developed during the TSC altimeter study are then used to establish the resulting 
design equation for measurement accuracy. As with the altimeter, measurement 
accuracy improves with S/N ratio up to a point, after which further increases 
in S/N provide little improvement (saturation effect! . In the wave spectro- 
meter, this effect occurs at unity signal-to-noise into the spectrum analyzer, 

(S/N) . Furthermore to achieve unity (S/N) . with a 2 dB modulation, it is 

isa 1 s a 

necessary to resolve to the order of a 19th of a (water) wavelength. 
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2*0 


HIGH RESOLUTION RADAR IMPLEMENTATION 


In addition to altimetry and wave height, a high resolution satellite 
radar may also be used to measure directional wave spectra. A nominal perfor- 
mance requirement for such an instrument is the measurement of spatial frequencies 
from 50 meters to 400 meters with an accuracy of 25%. This requirement is then 
compatible with a wave height measurement of 2.5 to 20 meters (peak-to-trough) 
with an accuracy of 25% since wave period is approximately twenty times peak-to- 
trough wave height. A directional resolution of at least (18°) is also assumed 
as a nominal requirement . 

There are at least three types of satellite radar Implementations 

capable in principle of achieving the above directional ocean wave spectra 

requirements. First, there is the high resolution side looking imaging radar. 

With this technique both image and directional wave spectra are continuously 

measured along the strip illuminated by the side looking antenna. The other two 

radar techniques provide a more coarse measurement of ocean wave spectra with a 

corresponding significant reduction in output data rate when compared to the 

radar imaging technique. The two frequency correlation technique of Ruck, 

Barrick and Kaliszewski ^ ^ and the short pulse wide-band radar technique dis- 

( 2 ) 

cussed by Tomiyasu are in fact mathematically equivalent. Both techniques 
perform a spectral analysis of the amplitude modulation (caused by ocean waves) 
on the return signal at a given analysis direction. A slow mechanically conical 
scanned antenna then permits the spectral measurement to be performed every 
A0° over the full 360° scan. 

The choice of implementing the short pulse wide-band radar technique 
is based primarily on the commonality of the implementation with that of the radar 
altimeter. In fact, a preliminary design analysis shows the resulting system 
parameters to correspond very closely to those of the altimeter design. The 
projected altimeter design can be fairly easily modified to perform ocean wave 
spectra measurements. The major modifications to the altimeter design which will 
allow it to measure directional surface spectra are: 
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1 4 The transmitter output is switched to the 2 meter antenna. 

The antenna beam is then steered 20° away from nadir and 
scanned conically with a period of 5.0 seconds- 

2. On receive the current full deramp processor is replaced 
by a surface wave pulse compressor. The compressed pulse 
output is detected and passed through a bank of 11 one- 
third (1/3) octave* filters which span the frequency range 
of 1 MHz to 9.0 MHz in logarithmic steps. The filter out- 
puts are then detected and integrated for 26 Misec per pulse. 

2 . 1 Basic Principles and Recommended Approach 

The basic purpose of this instrument is to measure the significant 
lower end (50-400 meters) of the gravity wave spectrum of the ocean surface as 
a function of angular direction. For an off nadir looking radar this requires 
the measurement of the low frequency portion of the spectrum of the mean power 
response of a matched filter receiver. This mean power response is the radar 
cross-section versus range, cr(R), which for a large enough off nadir look angle 
(20° for the present SEASAT-A geometry) will contain an amplitude modulation 
spectrum corresponding to the gravity wave spectrum. 

As described in section 2.0 of Part I of this report, a two frequency 
correlation technique measures the Fourier transform of a(R') as a function of 
frequency separation Af thus providing a direct measure of the desired amplitude 
modulation spectrum. A rigorous theoretical development of this technique for 
the measurement of ocean wave surface spectra is given in reference (1). With 
the short pulse radar technique (2), CT (R) is measured at the output of the 
matched filter receiver and hence must be further processed (spectrum analyzed) 
in order to obtain the desired amplitude modulation spectrum. As mentioned 
previously, both techniques are mathematically equivalent and simply represent 
different implementations of the same measurement process. The short pulse 
radar technique is chosen only because of the commonality of equipment with 


* A 1/3 octave filter has a bandwidth which is 25 % of its center frequency. 
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the short pulse altimeter. The following gives a rather simplistic view of 
the principles involved in the measurement of ocean surface spectra by means 
of a short pulse satellite radar. 


Consider a short pulse satellite radar looking at some angle 0 Q off 

nadir in the direction normal to the ocean wave (upwind or downwind) as illustrated 

in Figure 2.1. Then, as shown, there is a decrease in resolution along the 

C T 

surface Ax inversely proportional to the off nadir angle 0 q i.e. AR = and 

AX = - — f "" - n . Thus in order to faithfully produce an amplitude modulation on 
2 s in 0 

a(R) corresponding to the smallest spatial frequency of interest (50 meters) 
requires that AX ~ 25 meters. Note, however, that even if the pulse width T 
were a delta function; the amplitude modulation of a(R) can be distorted at 
large wave heights (peak-to-trough) if 0^ is made too small. For 20 meter peak- 
to-trough wave heights this consideration places a lower bound on 0 q at 8°. 

But the differential path length at the edge of the beam must also be on the 
order of the smallest wave period of interest. This length 


H (1 - cos 8 b/2 ) 

cos G sin 0 * 

o o 

which for a 2 meter antenna is about 50 meters at 0 = 20°. It is this con- 

o 

sideration then that places the lower bound on 0 q at about 20°. 

At this 20° off nadir angle the previous range resolution considera- 
tion (Ax ~ 25 meters) places a constraint on the maximum pulse length T at 
about 57 nsec. However, in order to have adequate S/N to detect small peak- 
to-peak amplitude modulations on the return it is required that 


B T - 
m c 


Y 


1 

2 


where the peak-to-peak modulation is given by 


±x 

1 - Y 


T c - compressed pulse width, 
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and is the bandwidth of the modulation. For ocean wave spectra the band- 
width is about equal to the frequency. Assuming the detection of a 2 dB 
modulation as a reasonable design goal, this consideration places the 
maximum pulse width at 6.0 nsec (for — 9.0 MHz corresponding to 50 meter 
wave periods). 

Under the conditions described above the output of the matched filter 
would appear as illustrated in Figure 2.2. The modulation of course would not 
have the regularity illustrated in the figure. Since spatial frequencies, &, 
of 50 to 400 meters are to be measured, the corresponding modulation frequencies 
of interest are from 9.0 MHz to 1.0 MHz since the relationship between the two 
is given by 

f = 

2 sin 9 l 
o 

Since the bandwidth of the modulation is approximately equal to f, a spectral 
analysis of this modulation by means of filters spaced at logarithmic intervals 
would seem to be a logical choice. With this selection, 11 filters having a 
1/3 octive bandwidth are sufficient to cover the 1.0 to 9.0 MHz frequency band 
of interest as shown in Figure 2.3. 

2 . 2 System Design Summary 

The preliminary design of a high resolution satellite radar for 
measuring directional ocean wave spectra is described in the following tables. 
This design meets all the specified performance requirements. These performance 
requirements are summarized in Table 2.1 and the system design is summarized in 
Table 2.2. The design presented in Table 2.2 was selected on the basis of 
compatibility with the altimeter design and as such probably does not represent 
an optimum configuration. A study to determine the trade offs involved in 
achieving an optimum configuration should be made during the initial phase of 
the program. A block diagram of the system is shown in Figure 2.4 and Table 2.3 
presents a preliminary power budget for this sensor. 
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Response (dB) 



Frequency (MHz) 


Figure 2.3 


Selected Filter Bank 


Table 2.1 


System Performance Requirements 


I. Ocean Wave Spectra 

Range 50-400 m 
Accuracy 25% 

II. Angular Directivity 

Range 360° 

Resolution 18° (5%) 

III. Detectable Modulation 2 dB 
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Table 2.2 
Design Summary 


I- 

Orbit Parameters 



a) Height 

725 km 


b) Inclination 

90° retrograde 


c) Eccentricity 

0.0064 maximum 

II. 

Radar Parameters 



a) Antenna beamwidth 

0.8° (2m dish) 


b) Pointing accuracy 

20 e = 0,1 


c) Antenna gain 

46.7 dB 


d) Peak power 

2 kW 


e) System losses (other than processing 
losses in pulse compressor) 

5 dB 


f) Noise figure 

5.5 dB 


g) Frequency 

13.9 GHz 


h) Uncompressed pulse width 

5 . 6 m>s 


i) Uncompressed pulse bandwidth 

180 MHz 


j) Compressed pulse width 

6.0 ns 


k) Compression ratio 

1000/1 


1) PRF 

> 400 Hz 


m) S/N (single pulse) 

3.0 dB 


n) Ocean cross section 

-10 dB 


o) Scan period 

5 sec 


p) Number of footprints in 360° 

154 

III. 

Linear FM Generation 



Type 

Surface wave 


Bandwidth 

180 MHz 


Pulse length 

5.6 (is 


Linearity of FM 

<0 . 27o 


Peak frequency deviation (one circle 
of variation across pulse) 

25 kHz 

IV. 

Pulse Compression 



Type 

Surface wave 


Bandwidth 

167 MHz 


Weighting 

-26 dB modified 
Taylor 


Processing loss 

0.55 dB 

V. 

Spectral Processing 



Type 

Filter band 


Number of filters 

11 


Frequency range 

1.0 to 9,0 MHz 


Filter spacing 

Logarithmic 


Filter bandwidth 

1/3 octave 


Integration time 

25 M>s 
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11 - 2*10 


PRF 

TRIGGER 


A/D kj 


2 KW PEAK POWER 



Figure 2,4 High Resolution Radar Block Diagram 




Table 2.3 

Signal-to-Noise Computation 



II-2-11 



3.0 


THEORETICAL ANALYSIS 


This section of the report describes the analysis of a directional 
wave spectrometer upon which much of the material presented in the previous 
section is based. While the overall approach and considerations presented 
here are believed to be correct, the following analysis should be considered 
a quick "first cut M at the problem because of the many simplifying and/or 
unsubstantiated assumptions. 

In section 3.1, the required receiver processing and underlying 
assumptions are established and used to develop design equations for the 
signal- to-noise ratio necessary for detection of the modulation caused by 
the gravity wave spectrum of the ocean surface. It is shown that the input 
to the spectrum analyzer consists of two noise processes, one being the 
desired modulation of 0 ° and the other due to the fact that the modulation 
is riding on a noise process whose mean is <j°. Standard statistical theory 
for the detection of a Gaussian signal in the presence of Gaussian noise is 
then used to develop an approximate expression for the resulting design 
constraint . 

In section 3.2, an expression for the standard deviation of the 

ocean wavelength measurement is derived. It is shown that the problem of 

measuring ocean wavelength is just the satellite altimeter (beam limited) 

problem in disguise. Results developed during the TSC altimeter study are 

then used to establish the resulting design equation for measurement accuracy. 

As with the altimeter, a saturation effect occurs when the "signal-to-noise" 

into the spectrum analyzer, (S/N) . , is near unity. Thus, a further increase 

isa 

in (S/N). beyond this value provides little improvement in measurement 

ISa 

accuracy. However, the ( s ^ N ^ iga is f° und to be a function of modulation 

index, Y, and the number of resolution cells per ocean wavelength. To achieve 

unity (S/N). with a 2 dB (peak-to-peak) modulation it is necessary to resolve 
i s a 

to the order of a 19th of a (water) wavelength. 
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IF 


3.1 Detectio n of Modulation on a Noise Process 

3.1.1 Assumed Processing 



B= 


At the output of the pulse compression filter the signal is given 


by: 


-2 

s(t) = [n + S + ySm(t) J [a (t)] cos(u>t + 0 (t)) 


(3.1) 


where 


N is the mean noise power 
S is the mean signal power 
m(t) is the modulation of the received signal power 
y is a modulation index 
a (t) is a Rayleigh process 
0 (t) is a random phase 


Then, a(t) cos 0 (t) and a(t) sin 0 (t) are independent Gaussian processes 
with bandwidth, B ( 
have unity power. 


with bandwidth, B c> equal to the pulse bandwidth. a(t) is assumed to 


Let B * be the bandwidth of the low pass filter. Assume that: 

-v 


But 


B C> V 


B m 


where B is the bandwidth of the modulation process m(t). 
m 


II-3-2 







After detection and low pass filtering: 


Z(t) = j S 2 (t-T)h(t)dt (3.2) 

» [N+S+YSm(t>] ^*a^(t-T)cos^((u(t-T)+0(t-i))h (t)dt 

where [N+S+YSm(t)] is taken outside the integral since it is assumed m(t) is 
constant over the response time of the filter (B c > B^) . 

The integral: 

q(t) = /a^(t-T)cos^(w(t-T) + 0 (t-T) )h(t)dt (3.3) 


Represents the result of passing the square of a narrow band Gaussian process 
(a(t)(cos u>t + 0 ( t ) ) through a low pass filter. The result of such an operation 
is to yield a Chi-square process with 2 (B^/B^) degrees of freedom. 

But the Chi-square process can be approximated by a normal process 
with mean 2(B /B p ) and variance 4(B /B p ) 

C *0 c *0 


Let 


where Tl(t) 
Thus 


B fiT 1 

q(t) « 2 ~ + 2 ^ Tl(t) 

■£- V b 

is a unit variance normal process. 


Z(t) « [(N+S) + YSm(t)] 


B 

2 — + 2 


= (N+S) 2 


ft) 


B 

+ 2YS ^ m(t) 



(3.4) 


(3.5) 

IKt) 


negligible 



That is, the cross product is negligible since 


Y < 1 


and 



Thus 


B 

/B \ 

/B? 


Z (t) » 2 (N+S) f + 2YS 

— |m(t) + 2 (N+S) 

f 

h 

TKt) 


(3.7) 


DC 


Signal 


Noise 


Thus using this approximation, the input to the spectrum analyzer 
consists of a DC bias plus two Gaussian noise processes (m(t) is simply assumed 
to be Gaussian). The process m(t) is the desired modulation, and the process 
T](t) is due to the fact that the modulation is riding on a (nearly white) noise 
process . 


Define B 

2YS — m(t) as "Signal" 


2 (N+S) 



T|(t) as "Noise" 


by: 


Then the signal to noise ratio into the spectrum analyzer is given 


(S/N) 



(3.8) 


Where S (to) 
m 

S^(to) is the spectral 


is the spectral density of the modulation process and 
density of the noise process. 
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Assuming band limited processes 


and 


Thus 


S(u>) 

tn 


B 

m 


(3.9) 


s n(»> 


1_ 

h 


(3.10) 


Recall that things are defined so that m and 11 have unit power. 



(3.11) 


or 



(3.12) 


3.1.2 Detection Criteria 

(3) 

Now the detection index for detecting Gaussian signal in the 
presence of Gaussian noise is: 

d = B T (S/N) 2 (3.13) 

a 

where B is the analysis bandwidth, and T is the total analysis time. i.e., 
a 

is the bandwidth of a filter in the spectrum analyzer, and T is the total 
integration time. Then 


where N 

P 




/ 2H9 sin0 \ 

b 2 ° 

' r» r A ' 


N 


o 

is the number of pulses averaged in the spectrum analyzer. 
(3) 

Good detection requires that 


(3.14) 


d > 25 
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Thus 


2H0 fi si 
2 

c cos 

represents a design constraint based on the systems capability to detect 
the modulation caused by ocean waves. 

3.2 Measurement of Ocean Wavelength 



n9 


N 


1,2 (nTs) (b 2 ) >5 

m 


(3.15) 


As was shown in the previous section, the input to the spectrum 
analyzer is a (nominally white) Gaussian process plus a process with a 
smaller bandwidth (Gaussian process) representing the modulation due to the 
wave structure, i.e,. 


Z(t) „ 2(N+S) 


(57) + 2yS (37) + 2 (N+S) V ^7 


nco 


In the spectral domain: 



the problem of measuring ocean wavelength is equivalent to locating the 
center of this hump in frequency. 
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Note this is just the satellite altimeter problem (beam limited) 
in disguise. In the spectral domain, what is required is to locate a hump 
of noise in the presence of uniform background noise. But this is just the 

(4) 

problem solved in the altimeter study, reference. (In fact, for a full 

stretch FM altimeter system, the altimeter problem is identical to this one 
since range becomes equivalent to frequency in a full stretch system.) 

This being the case, much of the following relies heavily on the 
results developed in the altimeter study including a consistent notation. 

Assume that a split-gate tracker is used to locate the center of 
the spectrum of the modulation as shown in Figure 3.1. Then the output of 
the tracker is given as 


e(w ) = S v - z v 

o eg i l g 1 


(3.16) 


where the are independent samples of the power, and the notation eg and 
■tg refer to early and late gate, respectively. 



Figure 3.1 Split-Gate Tracker 
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Further assume, for simplicity, that the shape of the modulation 
spectrum is symmetric about its center frequency. (The actual non- symmetric 
spectrum will of course cause a bias in this type estimate.) 


and 


Then 


E e(u> o ) = 0 


Now 


E e («> + Ato) 


= E v 

eg 1 


E V 

Ig i 


2Au> 


- B 


(S/N) 


isa 


a e(w o ) = 


V E V. V. + £ E V. V. _ 

T 1 J 1 J ^ 

*g eg 


(3.17) 


(3.18) 


= 2 


£ E V. (by symmetry) 

eg J 




(3.19) 


■ 2 { i 2 + ’i ’i} 


since E 


2 



for an exponential process. 


Hence 


a e(u) ) = 2 
v o' 




(3.20) 


II- 3-8 



Assuming the spectrum is nearly square: 


2 

i 


■ V 1 +S/N ls/ 


(3.21) 


2 V. = N E (1+S/N isa ) (3.22) 

i 


where N is the number of samples in the gate. 

Substituting eq. (3.21) and (3.22) into 3.20 yields 

cr 2 e(uu ) = 2(N_ + N_ 2 ) (1 + S/N. ) 2 
o h b isa 

Now the slope of the error signal (from eq. (3.18)) 

E e(u> o + Aw) 2 ( S/N ) lsa 

Aw - B Aw 

a 

Thus the error in estimating the center of the modulation hump is 


(3.23) 


(3.24) 


cr ~ 


OB 

e a 


(N e + N E 2 )B a 2 


» ~ 2 (S/N) 


1 + 


isa 


(S/N) 


(3.25) 


isa 


Now, from Figure 3.1, it is seen that the gates are assumed to cover the 
entire modulation spectrum, thus 


and 


N„B ~ B 
E a — m 


B B + B 
m a m 


w — 


1 + 


(S/N). 


isa J 


(3.26) 


(3.27) 
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Assuming that 


B « B 
a m 



B 

tn 




(3.28) 


Here again, the normalization of the error to the modulation bandwidth is 
somewhat analogous to that of the altimeter error which normalizes to wave 
height. Thus, for a given set of system parameters the error in estimating 
50 meter ocean wave lengths will be a factor of 8.0 greater than the error 
in estimating 400 meter ocean wave lengths. 


The error expression given by equation 
resolution basis. If more samples are averaged, 
decreases as ^/number. 


(3.28) is on a spectrum 
the standard deviation 


2H0 sinf) 

_ B N (see eq. 3.14) 

a 2 n p 

c cos 0 

o 


therefore 



(3.29) 


(3.30) 


Figure 3.2 shows the effect of (S/N)^ ga 
Note that a saturation effect occurs near (S/N)^ 3a 


on estimation accuracy. 
= 0 dB. 


1 1 - 3- 10 




(Const) 



<S/N* 8a 


S/N (dB) 


Figure 4.2 Estimation Accuracy vs (S/N). 

lScl 


Recall, however, that the "signal-to-noise" ratio is given by: 


S/fN isa Y ( S+N^ B 

m 


2 B, 


Thus, even for «> receiver signal-to-noise ratio, the effective 
signal-to-noise ratio becomes 

2 


o B 
.2 c 


B 


m 


B T 
m c 


(3.- 
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Hence, 

requires that 


to get 


in the neighborhood of unity 


"signal- to- noise " 

lSfl 


B T 
m c 



(3.32) 


This effect arises because you are trying to measure the modulation 

of 0°. But what you measure is a random return whose mean power is 0°. Thus 

you have to average several samples to estimate the mean power level. The 

less variation in a°, then the more accurately you must measure the mean 

power level to pick up small variations in it. But the only way to do this 

(unless you want to integrate every range resolution cell for several pulses) 

is to get several range cells packed into one wavelength on the water. (i.e. 

make B /B large), 
c m 

For Y = .23 (2 dB modulation) 


1 


B T 

m c 


> 19 


(3.33) 


Thus it is necessary to resolve to the order of a 19th of a (water) 
wavelength in order to recover the modulation reliably. 


Table 3.1 summarizes this requirement for various modulation indexes. 


Table 3-1 


Modulation (dB) 
(peak-to-peak) 

I 

(1/Y 2 ) 

2 

.23 

19 

4 

.43 

5.4 

6 

.60 

2.8 

8 

.73 

1.9 

10 

.83 

— 
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